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PREFACE 


This  report  is  the  third  of  a  series  dealing  with  the  manipulation 
and  interpretation  of  Landsat  digital  data.  The  concepts  and  procedures 
reported  herein  were  developed  under  the  In-House  Laboratory  Independent 
Research  Program,  Project  4A161101A91D,  Task  02,  Work  Unit  110-Q6,  "Use 
of  Automated  Procedures  for  Change  Detection  and  Land-Use  Mapping  from 
Satellite  Data,"  and  the  Earth  Resources  Satellite  Program,  Work  Unit 
31584,  "Programmed  Identification  and  Mapping  of  Forest  and  Habitat 
Categories." 

The  first  report  (TR  M-77-2,  Report  1)  of  this  series  discussed 
image  interpretation  and  the  theories  and  procedures  for  handling  re¬ 
flectance  geometry,  terrain  slope,  atmospheric  attenuation,  and  shadow 
effects.  These  extrinsic  effects  were  analyzed  and  an  analytic  expres¬ 
sion  for  their  correction  was  presented.  The  second  report  (TR  M-77-2, 
Report  2)  discussed  the  development  of  a  semiautomated  procedure  to 
classify  Landsat  radiance  data  in  terms  of  preselected  land-use  cate¬ 
gories  or  terrain  types. 

The  work  in  this  study  was  done  during  the  period  January  1977- 
September  1980  by  personnel  of  the  Environmental  Simulation  Branch  (ESB) 
under  the  direct  supervision  of  Mr.  H.  W.  West,  Project  Manager,  and 
Mr.  J.  K.  Stoll,  Chief,  ESB.  At  the  beginning  of  the  study,  the  ESB  was 
part  of  the  Environmental  Systems  Division  (ESD),  Mobility  and  Environ¬ 
mental  Systems  Laboratory  (MESL),  U.  S.  Army  Engineer  Waterways  Experi¬ 
ment  Station  (WES);  the  study  was  under  the  general  supervision  of 
Messrs.  B.  0.  Benn,  Chief,  ESD,  and  W.  G.  Shockley,  Chief,  MESL.  ESD 
later  became  part  of  the  Environmental  Laboratory  (EL)  of  WES;  the  study 
and  report  were  completed  under  the  general  supervision  of  Dr.  John 
Harrison,  Chief,  EL. 

Dr.  H.  Struve  (ESB)  developed  the  interactive  classification  pro¬ 
cedure  used  in  the  study.  The  report  was  prepared  in  FY  80  by  Ms.  M.  H. 
Smith  (ESB)  and  Dr.  Struve.  Mr.  A.  N.  Williamson,  Computations  and 
Analysis  Group,  Mobility  Systems  Division,  Geotechnical  Laboratory, 


produced  the  graphics  material  of  the  Landsat  digital  data  used  in  the 
preparation  of  this  report. 

COL  John  L.  Cannon,  CE,  and  COL  Nelson  P.  Conover,  CE,  were  Com¬ 
manders  and  Directors  of  WES  during  the  period  of  study  and  report 
preparation.  Mr.  F.  R.  Brown  was  Technical  Director  during  this  period. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con¬ 
verted  to  metric  (SI)  units  as  follows: 


Multip ly 


To  Obtain 


feet 

inches 

miles  (U.  S.  customary) 


0.3048 

25.4 

1.609344 


metres 

millimetres 

kilometres 
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ACQUISITION  OF  TERRAIN  INFORMATION  USING 
LANDSAT  MULTISPECTRAL  DATA 


APPLICATION  OF  AN  INTERACTIVE  CLASSIFICATION 
PROCEDURE  IN  SOUTH  LOUISIANA 

PART  I:  INTRODUCTION 

Background 

1.  The  capability  of  a  Landsat  satellite  to  view  the  earth's  sur¬ 
face  from  a  vantage  point  900  to  950  km  above  the  earth  has  revolution¬ 
ized  the  concept  of  data  collection  for  both  military  and  Civil  Works 
purposes.  Not  only  does  the  "eye  in  the  sky"  see  from  this  perspective, 
but  it  consistently  follows  its  track,  circling  the  earth  with  such  pre¬ 
cision  that,  at  regular  intervals,  it  views  (in  a  continuum)  approxi¬ 
mately  the  same  strip  of  the  earth's  surface.  On  board  the  satellites 
of  the  Landsat  Observation  System,  multispectral  scanners  (MSS)  and  re¬ 
turn  beam  vidicon  cameras  (RBV)  continually  sample  the  spectral  reflec¬ 
tances  from  the  earth's  surface  and  record  them  on  magnetic  tape  for 
analysis  and  image  production.  These  tapes  also  contain  ephemeris  data 
and  other  information  about  the  environment  at  the  time  of  collection 
that  permit  reasonable  assessment  of  the  quality  of  the  data  from  the 
earth's  surface.  Thus,  the  Landsat  System  puts  a  mass  of  current,  reli¬ 
able  data,  collected  by  consistent  procedures,  in  the  hands  of  geogra¬ 
phers,  strategists,  economists,  and  researchers  for  use  in  their  search 
for  expanded  knowledge  of  the  surface  of  the  earth,  both  in  our  nation 
and  around  the  world. 

2.  On  the  domestic  scene,  in  addition  to  the  obvious  applications 
of  Landsat  data  analysis  for  defense  strategy,  analyses  of  the  data  hold 
promise  for  monitoring  the  use  of  the  nation's  waterways  and  coastal  re¬ 
gions,  monitoring  waste  disposal,  assessing  flood  damage,  locating  urban 
sprawl  and  industrial  development,  and  monitoring  change  in  land  use. 

3.  Potential  worldwide  military  applications  are  for  military 


7 


support  information  and  for  building  data  banks  to  be  used  in  future  re¬ 
search  or  in  times  of  emergency.  Nonmilitary  needs  that  could  be  served 
by  Landsat  data  analyses  are  in  the  area  of  protection  of  water  re¬ 
sources  and  monitoring  the  changes  made  by  the  exploitation  of  other 
natural  resources.  The  emerging  third  world  nations  struggling  to 
become  responsible  and  self-sufficient  could  be  served  by  knowledge  of 
their  nations'  land,  its  use,  and  changes  brought  about  by  various  land- 
management  methods. 

4.  Landsat  data  have  been  available  only  a  very  short  time,  the 
first  satellite  being  launched  in  July  1972.  There  are  now  three  satel¬ 
lites,  the  second,  Landsat  2,  launched  in  January  1975,  and  the  third, 
Landsat  3,  launched  in  early  1978.  The  first  two  sampled  reflectance 
values  from  the  earth's  surface  by  six  sets  of  four  sensors  each  in  the 
spectral  range  of  0.5  through  1.1  pm.  The  third  satellite  carries,  in 
addition  to  these  four  sensor  sets,  a  fifth  sensor  in  two  of  the  six  sets 
that  operates  in  the  thermal  infrared  region  from  10.4  to  12.6  pm.  This 
sensor,  operating  in  the  fifth  spectral  reflectance  band,  will  result  in 
one-third  the  resolution  of  the  other  four  spectral  bands.  Computer 
compatible  tapes  (CCT's)  obtained  from  Landsat  3  will  not  carry  informa¬ 
tion  from  the  new  band  five  until  the  data  have  been  thoroughly  evalu¬ 
ated.  Early  indications  suggest  possible  distortion  of  data  beyond  re¬ 
covery,  preventing  the  availability  of  these  kinds  of  data  until  the 
launching  of  another  satellite.  Additional  information  on  the  operation 
and  products  of  satellites  in  the  Landsat  Observation  System  is  available 
in  published  reports.* 

5.  Many  investigators,  private  as  well  as  government  agencies, 
have  been  occupied  with  the  analysis  of  Landsat  data.  Among  the  govern¬ 
ment  agencies,  the  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES) 


*  National  Aeronautics  and  Space  Administration.  1976.  "Landsat  Data 
Users  Handbook,"  Document  No.  76SD4258,  Goddard  Space  Flight  Center, 
Greenbelt,  Md . 

Jaffe,  L.  1975.  "Inputs  Requested  from  Earth  Resources  Remote  Sens¬ 
ing  Data  Users  Regarding  Landsat-C  Mission  Requirements  and  Data 
Needs,"  Attachments  A-E,  National  Aeronautics  and  Space  Administra¬ 
tion,  AN-OA-76-B,  Goddard  Space  Flight  Center,  Greenbelt,  Md. 
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has  taken  keen  interest  in  the  potential  of  Landsat  data  to  identify 
water  bodies,  land  use,  and  wetlands;  and  to  detect  the  transport  and 
deposition  of  suspended  material  in  streams  and  rivers,  the  illegal  use 
of  marginal  river  and  coastal  boundaries,  and  the  changes  in  land  use 
caused  by  man  and  by  natural  phenomena.  In  these  endeavors,  the  analysis 
of  Landsat  data  indicates  its  potential  for  monitoring  land  use  and  its 
changes.  However,  the  analysis  procedures  are  still  in  an  almost  embry¬ 
onic  state.  The  information  resides  in  the  analogs,  and  the  digital 
transforms  are  stored  on  CCT's  for  users  to  interpret  in  their  investi¬ 
gative  studies.  While  much  has  been  accomplished  by  National  Aeronautics 
and  Space  Administration  (NASA)  researchers  in  their  transfer  of  data  to 
the  CCT's,  interpretation  remains  for  the  accomplishment  of  the  user. 

The  first  two  reports*  in  the  "Acquisition  of  Terrain  Information  Using 
Landsat  Multispectral  Data"  series  describe  a  procedure  for  transforming 
and  interpreting  Landsat  data. 

6.  The  first  report  (TR  M-77-2,  Report  1)  discussed  image  inter¬ 
pretation  and  the  theories  and  procedures  for  handling  reflectance 
geometry,  terrain  slope,  and  atmospheric  and  shadow  effects.  These  ex¬ 
trinsic  effects  were  analyzed  and  a  consolidation  equation  for  their  cor¬ 
rection  was  presented.  The  equation  can  be  used  to  compute  terrain  radi¬ 
ance  as  a  function  of  wavelength  as  measured  by  the  Landsat  Sensor  System 
in  watts  per  square  centimetre  per  micrometre  (Wcm  pm  ).  Plates  1-4, 
taken  from  Struve,  Grabau,  and  West  (1977a),  show  a  generalized  diagram 
of  procedural  steps  to  interpret  a  Landsat  image  and  explain  the  de¬ 
tailed  steps  in  atmospheric,  slope,  and  shadow  correction  procedures. 

7.  The  second  report  (TR  M-77-1,  Report  2)  discussed  the 

*  Struve,  H. ,  Grabau,  W.  E.,  and  West,  H.  W.  1977a.  "Aquisition  of 
Terrain  Information  Using  Landsat  Multispectral  Data;  Correction  of 
Landsat  Spectral  Data  for  Extrinsic  Effects,"  Technical  Report  M-77-2, 
Report  1,  U.  S.  Army  Engineer  Waterways  Experiment  Station,  CE,  Vicks¬ 
burg,  Miss. 

Struve,  H. ,  Grabau,  W.  E.,  and  West,  H.  W.  1977b.  "Aquisition  of 
Terrain  Information  Using  Landsat  Multispectral  Data;  An  Interactive 
Procedure  for  Classifying  Terrain  Types  by  Spectral  Characteristics," 
Technical  Report  M-77-2,  Report  2,  U.  S.  Army  Engineer  Waterways  Ex¬ 
periment  Station,  CE,  Vicksburg,  Miss. 


development  of  a  semiautomated  procedure  to  classify  Landsat  radiance 
data  in  terms  of  preselected  land-use  categories  or  terrain  types.  The 
procedure  can  be  applied  directly  to  the  primary  data  on  the  CCT's,  as¬ 
suming  that  the  relations  among  the  primary  data  are  the  same  as  those 
among  processed  data;  i.e.,  classification  results  would  not  be  improved 
if  applied  to  processed  or  corrected  data  instead  of  primary  data. 

Actual  data  from  the  Mississippi  River  floodplain  north  of  Vicksburg 
(Satartia  site)  were  used  in  the  development  process.  The  interactive 
procedure  described  in  Struve,  Grabau,  and  West  (1977b)*  was  very  suc¬ 
cessful  in  the  classification  of  water  and  forest  pixels  in  the  data 
from  the  demonstration  site.  Computer  programs  developed  in  this  pro¬ 
cedure  can  be  used  easily  by  engineers  or  trained  technicians  in  an  in¬ 
teractive  fashion  to  classify  water  and  forest  from  Landsat  data,  pro¬ 
vided  some  preselected  guidance  is  available.  Plates  5  and  6,  taken 
from  Report  2,  describe  the  procedural  steps  for  a  generalized  guided 
classification  of  Landsat  data  and  detailed  steps  in  a  procedure  for  the 
selection  of  classification  criteria. 


Objective  and  Scope 


8.  The  overall  objective  of  the  study  series  is  to  develop  and 
demonstrate  automated  procedures  for  obtaining  land-use  classifications 
from  Landsat  data  that  are  acceptable  for  use  by  the  Corps  Divisions  and 
Districts  in  detecting  land-use  change  and  producing  land-use  maps.  The 
specific  objective  of  this  study  was  to  demonstrate  the  application  of 
the  guided  classification  procedure  (in  locations  other  than  the  Satartia 
site)  to  identify  water,  forested,  and  nonforested  land  use;  to  explore 
the  possibility  of  extending  the  procedure  to  include  urban  and  open 
space  classes;  and  to  investigate  an  unguided  classification  procedure. 


9.  Landsat  data  from  two  flights  over  Baton  Rouge,  La.,  and 


*  Ibid. 
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vicinity  were  studied:  an  older  flight  in  October  1972  (tapes  on  hand 
at  WES)  and  a  more  recent  flight  in  February  1977  (tapes  ordered  from 
the  Earth  Resources  Observation  System  (EROS)  Data  Center*).  February 
1972  aerial  photography  was  used  with  the  October  1972  overflight,  and 
February  1974  aerial  photography  was  obtained  for  the  February  1977  over¬ 
flight.  At  the  time  of  the  study,  1974  photography  was  the  latest  avail¬ 
able.  The  procedural  steps  outlined  in  Plates  5  and  6  were  followed  to 
identify  water,  forest,  and  other  (nonforest,  nonwoods)  land  use  in 
the  Baton  Rouge  vicinity.  Automated  procedures  were  developed  to  trans¬ 
late  the  data  from  satellite  magnetic  tape  to  a  standard  WES  format  for 
interpretation,  to  extract  data  for  specific  geographical  locations,  and 
to  sort  data  into  unique  signatures  and  count  their  occurrence.  Pro¬ 
cedures  to  rectify  Landsat  data  to  aerial  photography  and  an  automated 
procedure  to  rectify  Landsat  data  to  Universal  Transverse  Mercator  (UTM) 
coordinates  were  demonstrated.  Both  guided  and  unguided  procedures  were 
used  to  classify  the  Landsat  data.  A  comparison  of  the  land-use  classes 
defined  by  the  Landsat  procedures  and  aerial  photography  interpretation 
was  made  for  water,  forest  (or  woods),  urban,  and  open  space  (e.g., 
fields).  Also,  the  land  use  interpreted  from  1977  Landsat  data  was  com¬ 
pared  with  1972  Landsat  data. 


*  U.  S.  Department  of  the  Interior,  Geological  Survey,  EROS  Data 
Center,  Sioux  Falls,  S.  Dak. 


PART  II:  DESCRIPTION  OF  THE  STUDY 


Study  Sites 


Baton  Rouge  and  vicinity 

10.  Baton  Rouge,  La.,  and  vicinity  (Figure  1)  was  selected  for 
study  for  the  following  reasons. 

<1.  It  contains  large  areas  of  the  land-use  classes  of 
interest  in  this  study. 

It  contains  distinctive  landmarks  such  as  highway  inter¬ 
sections,  water  bodies,  and  airport  runways  for  use  as 
reference  points. 

£.  There  are  several  Landsat  cloud-free  scenes  containing 
the  region. 

<1.  Two  universities  and  numerous  State  and  Federal  govern¬ 
ment  offices  are  available  for  information. 

It  is  less  than  150  miles*  from  WES. 

11.  The  primary  reason  for  selection  was  the  presence  of  large 

areas  of  water,  forested  (or  wooded),  urban,  and  open  areas  (Figure  2). 

The  Mississippi  River  occupies  the  center  of  the  region,  top  to  bottom, 

with  several  lakes  and  a  portion  of  the  Intracoastal  Waterway  included. 
Baton  Rouge  is  on  the  east  side  of  the  river  and  contains  a  variety  of 

land  uses  found  typically  in  urban  areas.  There  are  oil  refineries  and 

other  industrial  areas,  business  sections,  new  and  old  residential  areas 
college  campuses,  and  parks.  On  the  west  side  of  the  river  is  an  iso¬ 
lated  oil  refinery  and  a  smaller  urban  area  (Port  Allen)  consisting 
almost  entirely  of  residences.  A  large  wooded  area  surrounds  the  Intra¬ 
coastal  Waterway  in  the  southwest  portion  of  the  region.  The  remainder 
of  the  region  west  of  the  river  is  open  space  covered  by  long,  narrow 
sugarcane  fields,  scattered  aggregates  of  dwellings,  crossroad  stores, 
and  a  small  oil  field  or  two. 


*  A  table  of  factors  for  converting  U.  S.  customary  units  of  measure¬ 
ments  to  metric  (SI)  units  is  presented  on  page  6. 
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Location  of  study  site 


Figure  2.  Baton  Rouge  and  vicinity,  February  1972  aerial  photography 


Port  Allen 


12.  Port  Allen,  a  subset  of  the  Baton  Rouge  site  including  an  oil 
refinery,  woods,  and  fields,  was  chosen  first  for  study  because  of  the 
variety  of  land  uses  in  a  smaller  area  and  because  the  concentration  of 
houses  in  Port  Allen  offered  an  opportunity  for  better  understanding  of 
reflectance  from  residential  areas.  February  1972  aerial  photography 
was  used  for  land-use  sample  locations  in  interpretation  of  October  1972 
Landsat  data. 

13.  The  classification  for  the  Port  Allen  area  was  used  in  expand¬ 
ing  the  study  area  to  the  Baton  Rouge  site.  The  larger  site  was  also 
classified  using  1977  Landsat  data  with  land-use  interpretation  from 
1974  aerial  photography. 

Atchafalaya  Basin 

14.  A  site  in  the  Atchafalaya  Basin  was  used  as  a  demonstration  of 
the  automated  procedure  to  rectify  Landsat  data  to  UTM  coordinates  and  to 
compare  the  success  of  the  classification  procedure  to  a  previous  manual 
land-use  classification  by  interpretation  of  aerial  photography. 

Data  and  Equipment 

15.  Landsat  data  from  two  overpasses  were  used  in  this  study,  one 
on  1  October  1972  and  one  on  28  February  1977.  Information  about  each 
overpass  is  included  in  Table  1.  The  images  produced  from  data  in  band  5 
of  each  overpass  are  shown  in  Figures  3  and  4.  The  1972  data  were  col¬ 
lected  by  Satellite  1  and  the  1977  data  by  Satellite  2.  The  CCT's  were 
obtained  from  the  EROS  Data  Center.  The  CCT's  and  their  formats  were 
discussed  in  Report  1  of  this  study  series.*  For  details  of  Landsat 
data  records,  see  the  users  handbook.** 

16.  The  1972  aerial  photography  was  black  and  white  photography 
acquired  from  the  U.  S.  Department  of  Agriculture  flown  at  20,000  ft 
with  a  6-in.  focal  length  camera.  An  enlargement,  1:12,000,  was  used 
for  land-use  interpretation.  The  1974  aerial  photography  was  color  in¬ 
frared  photography  acquired  from  EROS  Data  Center  flown  by  NASA  aircraft 

*  Struve,  Grabau,  and  West  (1977a),  op.  cit. 

**  NASA,  op.  cit. 
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at  131,500  ft.  An  enlargement,  1:65,500,  was  used  in  land-use 
interpretation . 

17.  A  Honeywell  G-635  computer  was  used  to  process  the  data.  The 
computer  consists  of  a  256  K  word  core  memory,  disc  and  magnetic  tape 
storage  subsystems,  two  line  printers,  and  three  Data  Net  30  processors. 
The  normal  processing  of  Landsat  data  is  similar  to  other  array  pro¬ 
cessing  techniques  by  sequential  procedures.  One  or  more  records  (scan 
lines  of  data)  are  read  sequentially  into  a  buffer  from  a  data  tape. 
Sequential  processing  of  each  element  in  the  record  takes  place  with  the 
resulting  elements  placed  in  an  output  buffer.  When  the  buffer  is 
filled,  the  data  are  written  sequentially  to  an  output  tape  or  to  a 
printer  for  hard-copy  records.  The  input  buffer  is  refilled  and  the 
process  continues  until  all  data  are  processed. 

18.  A  line  printer  is  used  to  print  digital  values  extracted  from 
the  Landsat  CCT.  A  program  has  been  developed  to  print  the  extracted 
data  in  swaths  that,  when  joined,  present  a  digital  picture  of  the  de¬ 
sired  portions  of  the  Landsat  scene. 

19.  A  film  reader/writer  (described  in  Appendix  A  of  Struve, 
Grabau,  and  West  (1977b)*)  is  used  to  present  a  picture  in  shades  of 
gray  to  represent  the  data.  A  light  beam,  with  intensity  varying  as  the 
input  data  values  vary,  moves  across  a  piece  of  film  mounted  on  a  rotat¬ 
ing  drum.  The  light  exposes  a  minute  square  of  film  producing  a  shade 
of  gray  determined  by  light  intensity.  Every  data  value  is  represented 
by  a  small  gray  square.  When  the  squares  are  all  produced,  a  "picture" 
of  the  data  can  be  seen  in  shades  of  gray. 

20.  Since  the  film  writer  is  restricted  to  the  use  of  square 
elements  in  the  picture,  a  procedure  was  developed  to  stretch  the  pic¬ 
ture  along  the  flight  line  of  the  satellite  to  more  nearly  represent 
the  actual  earth  surface  viewed  by  the  sensor.  This  was  accomplished 
for  the  site  in  this  study  by  repeating  every  third  and  twentieth  scan 
line  of  data.  Sites  in  latitudes  significantly  different  would  require 

a  different  stretch  or  some  geometric  pixel  size  normalization  procedure. 


*  Struve,  Grabau,  and  West  (1977b),  op.  cit. 
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21.  The  film  writer  is  used  to  produce  pictures  of  the  land-use 
classes  as  interpreted  from  the  Landsat  data.  The  input  to  the  film 
writer  is  a  magnetic  tape  produced  from  the  Landsat  data  array  with 
class  codes  replacing  the  digital  reflectance  values  of  each  pixel. 

This  tape  is  read  into  the  film  writer,  which  produces  a  unique  shade  of 
gray  for  each  land-use  class  resulting  in  a  map  of  land  use.  The  film 
writer  products  are  black  and  white  pictures  for  study,  halftones  for 
report  or  display  products,  or  color  masks  for  color  prints. 

Classification  Procedures  of  the  Subsite  with  1972  Data 


Landsat  data 

22.  The  procedures  outlined  in  Plates  5  and  6  were  followed  to 
interpret  and  classify  the  Landsat  data.  Steps  1  and  2  were  reversed  in 
order  to  check  the  coverage  and  quality  of  the  data  before  incurring  the 
expense  of  acquiring  CCT's  (normally  $200  for  a  scene).  The  images  of 
the  two  Landsat  overpasses  were  obtained.  The  location  and  quality  were 
verified  (Table  1),  and  the  CCT’s  were  ordered  from  EROS  Data  Center. 

The  area  of  interest  was  located  on  the  Landsat  image,  and,  with  the  use 
of  a  Gerber  variable  scale,  the  corners  were  identified  in  terms  of 
pixels  and  scan  lines.  Figure  5  shows  the  orientation  of  the  Landsat 
data  for  the  site  on  a  topographic  map,  scale  1:250,000.  Since  the 
satellite  travels  at  an  angle  of  9  to  11  deg  from  north  in  the  latitude 
of  Baton  Rouge,  the  Landsat  data  including  the  area  of  interest  are 
tilted  at  the  Landsat  orbit  angle  with  the  scan  lines  perpendicular  to 
the  line  of  flight. 

23.  The  CCT's  were  read  and  unpacked  from  the  Landsat  format  into 
a  more  convenient  format  for  use  on  the  36  bit  word  Honeywell  G-635  com¬ 
puter.  This  format,  called  WES  format,  places  certain  header  informa¬ 
tion  on  the  file  and  places  data  from  one  pixel  into  one  computer  word. 

A  pixel  is  understood  to  be  the  spot  of  earth  (approximately  57  by  79  m) 
viewed  by  a  Landsat  sensor  at  exposure  time  or  the  numerical  value  repre¬ 
senting  the  reflectance  from  that  earth  spot.  All  the  data  from  the  west 
half  of  a  scene  reside  on  one  WES-formatted  tape,  that  from  the  east  half 
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Figure  5.  Topographic  map  showing  UTM  coordinates  and  orientation 
of  extracted  1  October  1972  Landsat  data  containing  area  of 

interest 
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on  another.  The  data  from  the  area  of  interest  were  extracted  from  the 
WES-formatted  tapes  and  were  written  on  a  new  tape  in  the  same  format. 

24.  For  a  visual  concept  of  the  data,  pictures  of  the  site  as 
sensed  by  each  band  were  produced  with  the  film  writer  by  grouping  the 
reflectance  values  into  nine  groups  (Figure  6).  Each  group  was  then 
assigned  a  distinct  shade  of  gray.  Low  reflectance  values  were  assigned 
light  shades  and  high  values  were  assigned  dark  shades.  The  Mississippi 
River  and  lakes  appear  very  light  in  band  7  compared  to  other  surfaces 
due  to  their  absorption  of  spectral  wavelengths  0.8  to  1.1  pm.  Bands  4 
and  5  show  greater  contrasts  of  nonwater  surfaces  than  bands  6  and  7, 
which  should  make  them  more  useful  for  classifying  forest  and  other  non¬ 
water  surfaces.  The  lightest  areas  in  bands  4  and  5  are  forested  areas. 
The  flat-appearing  area  (next  lightest)  in  band  4  to  the  west  of  the 
river  in  the  center  portion  of  the  scene  is  occupied  mostly  with  sugar¬ 
cane  fields  and  other  open  space.  This  area  in  the  other  bands  appears 
to  contain  a  wider  range  of  values,  which  would  increase  the  difficulty 
of  separating  open  spaces  and  urban  areas.  There  also  appears  to  be  a 
change  in  ratio  of  one  land-use  class  to  another  from  band  to  band;  e.g. 
the  ratio  of  forest  to  water  appears  to  be  quite  different  in  band  4 
from  that  in  band  7. 

Geometric  accordance  of  Landsat 
data  with  aerial  photography 

25.  The  guided  classification  procedure  requires  some  knowledge  of 
the  land-use  classes  as  close  to  the  time  of  the  satellite  overpass  as 
possible.  Samples  are  located  on  the  available  aerial  photography  and 
reflectance  values  from  each  band  are  extracted  from  the  corresponding 
locations  in  the  Landsat  data.  To  locate  the  correct  sample,  the  Landsat 
data  must  be  brought  into  geometric  accordance  with  the  aerial  photogra¬ 
phy.  This  was  done  by  locating  three  points  that  were  easily  identified 
in  both  data  sets  to  be  used  as  transfer  points.*  For  the  Port  Allen 
site  (Figure  7),  the  three  points  A,  B,  and  C  were  easily  identified  in 
the  band  5  Landsat  data  (Figure  8).  Points  A  and  B  are  intersections  of 
clearly  defined  highways,  and  point  C  is  identified  by  the  high 


Figure  6.  Gray-scale  pictures  from  Landsat  scene,  28  February  1977, 
using  Landsat  reflectance  values  separated  by  bands 
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Figure  8.  Transfer  point  locations  for  geometric  accordance  of  1972 
Landsat  with  1972  aerial  photography  data 

reflectance  from  the  northwest  corner  tank  in  rows  of  large  oil  storage 
tanks.  The  origin  was  assumed  at  A,  and  the  distances  of  B  and  C  from  A 
with  the  pixel  and  scan  numbers  of  the  corresponding  points  in  the  Land- 
sat  data  supplied  the  necessary  information  for  the  solution  of  the  fol¬ 
lowing  equations: 
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0£  =  Landsat  flight-path  angle  with 
respect  to  true  north 

=  number  of  X-coordinate  pixels  of 
the  transfer  point  within  the 
Landsat  data  set 

Dv.  =  width  of  the  Landsat  pixel 

A 1 

Ny .  =  number  of  Y-coordinate  pixels  of 
the  transfer  point  within  the 
Landsat  data  set 

Dy .  =  length  of  the  Landsat  pixel 

Equations  1  and  2,  representing  each  respective  trans¬ 
fer  point,  were  written  with  the  appropriate  substi¬ 
tutions  and  then  simultaneously  solved  for  the  flight 
path  angle  0  and  the  size  (D^  and  Dy)  of  the  pixel 

within  each  Landsat  data  set.  The  average  values  for 

0  ,  D  ,  and  D  were  then  computed  and  used  in  the 

I  A  i 

equations  above  to  bring  the  terrain  reference  samples 
and  CCT  value  data  into  geometric  accordance."* 

26.  These  transfer  constants  were  then  used  to  transfer  the  loca¬ 
tion  of  land-use  samples  to  the  Landsat  pixels  representing  the  spectral 
value  of  those  samples.  A  listing  of  the  digital  Landsat  band  5  data 
from  the  site  area  was  obtained.  For  ease  of  selection,  the  Landsat  data 
had  been  stretched  (repeating  third  and  twentieth  scan  lines)  to  present 
the  data  in  a  more  easily  recognized  relationship  to  the  photography. 

This  stretching  by  inserting  pixels  over  the  same  distance  shortens  the 
pixel  length  producing  square  pixels,  =  Dy  .  Table  2  lists  the  ac¬ 

cordance  data  for  1  October  1972  Landsat  data  to  February  1972  aerial 
photography  and  28  February  1977  Landsat  data  to  February  1974  aerial 
photography  for  Baton  Rouge  and  vicinity.  A  computer  program  to  deter¬ 
mine  these  constants  using  unstretched  data  was  developed  and  later  in¬ 
serted  into  the  procedure.  The  program  requires  at  least  five  transfer 
points  over  an  area  including  the  study  area.  This  automated  program 
was  later  used  for  rectifying  classified  portions  of  the  Landsat  scene 
to  UTM  coordinates  on  topographic  maps.  Table  3  lists  the  data  for  rec¬ 
tifying  sites  to  UTM  coordinates  on  a  1:250,000  topographic  map. 

*  Ibid. ,  pp  53-54. 
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Land-use  sample 
locations  and  data 


27.  Twenty-one  land-use  samples  were  selected  from  the  February 
1972  aerial  photograph  (Figure  7).  Three  samples  (1-3)  were  selected 
from  wooded  areas;  three  were  in  water  (4-6);  six  were  in  the  urban 
residential  area,  Port  Allen  (7-12);  three  were  in  an  urban  industrial 
area  (13-15);  and  six  (16-21)  were  in  open  space.  Equations  1  and  2 
were  used  with  the  appropriate  rectification  constants  and  distances  of 
the  sample  locations  from  transfer  point  A  to  determine  the  correspond¬ 
ing  location  (pixel  and  scan)  in  the  Landsat  data. 

28.  Reflectance  values  of  each  terrain  sample  and  their  eight 
surrounding  terrain  pixels  were  located  in  a  digital  listing  of  the 
Landsat  data  array.  An  example  of  the  3  by  3  arrays  of  data,  Sample  1 
(woods),  is  listed  below. 


Band  4 

22  22  22 
20  21  24 
20  22  22 


Band  5 

13  13  12 
13  11  13 
12  12  12 


Band  6 

34  32  32 

33  33  35 

34  34  34 


Band  7 

21  21  21 
22  22  21 
22  20  22 


Spectral  data  analysis 
of  Landsat  1972  data 

29.  The  mean  and  standard  deviation  were  computed  for  each  3  by 
3  data  array  to  determine  if  the  samples  were  representative  of  the  ter¬ 
rain  type  from  which  they  were  selected  on  the  aerial  photo.  The  three 
samples  from  woods  (samples  1-3  in  Table  4)  had  standard  deviations  less 
than  1.5.  Struve,  Grabau,  and  West  (1977b)*  established  a  standard  devi¬ 
ation  of  1.5  or  less  as  a  fair  criterion  for  homogeneity  of  a  sample; 
therefore,  these  samples  were  retained  in  the  classification  procedure. 
The  water  samples,  samples  4-6  (Table  4),  also  contained  none  with  devi¬ 
ations  higher  than  1.5.  Most  of  the  other  3  by  3  arrays  have  standard 
deviations  much  higher  than  the  homogeneity  criterion.  This  was  to  be 
expected  from  the  nature  of  these  kinds  of  land  use.  The  mixtures  of 


*  Struve,  Grabau,  and  West  (1977b),  op.  cit. 
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structures,  trees,  and  open  spaces  in  urban  areas  and  the  variety  of 
planted  fields,  open  pastures,  wet  areas,  and  plowed  fields  in  open 
spaces  present  a  "checkerboard"  of  high-  and  low-valued  pixels  in  both 
urban  areas  and  in  open  spaces.  With  the  exception  of  the  older  residen¬ 
tial  areas,  the  band  4  reflectance  from  urban  areas  is  generally  higher 
than  that  from  open  space  (Figures  2  and  6).  For  this  reason,  it  was 
believed  that  urban  could  be  separated  from  open  space  using  band  4,  and 
the  samples  were  kept  in  the  analysis  even  though  they  did  not  meet  the 
homogeneity  criterion. 

30.  Using  the  means  of  the  four  bands  as  representative  signatures 
of  the  sample  types,  a  similarity  analysis  was  made  with  the  clustering 
technique  used  by  Struve,  Grabau,  and  West  (1977b).*  The  signatures 
with  an  allowable  deviation  were  entered  into  the  computer  program. 
Signatures  were  then  grouped  by  those  falling  within  the  allowable  range 
of  each  band.  The  range  of  allowable  deviation  was  increased  until 
unique  clusters  were  formed,  i.e.,  samples  did  not  overlap  clusters.  The 
results  of  this  interactive  cluster  technique  for  Port  Allen  La.  isat  1972 
data  are  shown  in  Figure  9.  A  black  dot  indicates  similarity  of  the  two 
samples;  e.g.,  sample  1  is  similar  to  samples  2  and  3.  Samples  20  and 
21  were  found  to  be  more  similar  to  urban  (new  residential)  than  to  open 
space,  so  they  were  eliminated  from  the  classification  procedure.  Be¬ 
cause  it  did  not  encroach  upon  another  class,  sample  12  was  kept  with 
the  urban  (residential)  samples  even  though  it  was  not  similar.  The 
higher  industrial  urban  reflectance  samples  clustered  separately  from 
the  residential  urban  samples.  Their  range  is  shown  as  a  double  bar  on 
the  bar  plots  of  terrain  classes  in  Figure  10.  In  Figure  10,  plots  a- 
d  are  plots  of  the  reflectance  ranges  of  the  clusters  in  Figure  9  by 
bands  4-7.  These  cluster  ranges  were  derived  from  the  cluster  analysis 
with  a  deviation  of  1.75  and  are  listed  in  Table  5.  The  range  (mean 
+1.75)  is  listed  for  each  band  of  each  sample  in  the  cluster.  The 
cluster  range  was  then  determined  as  the  range  from  lowest  to  highest 
reflectance  in  the  cluster.  From  the  bar  plots  of  these  ranges 


Figure  9.  Cluster  analysis  of  1  October  1972  Landsat  data 
from  Port  Allen,  1.75  variance 


Figure  10),  it  can  be  seen  that  water  is  separated  completely  from  other 
clusters  in  band  7;  woods  are  separated  in  bands  4  and  5.  There  is  no 
clear  separation  of  urban  clusters  from  open  space  anywhere  except  for 
those  clusters  with  high  reflectance  in  bands  4  and  5.  Since  a  greater 
proportion  is  separated  in  band  4,  it  was  used  to  establish  the  urban 
class  limits.  There  is  no  chance,  with  this  procedure,  to  separate  lower 
reflectance  in  urban  areas  from  reflectance  in  open  space. 
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31.  The  following  ranges  were  accepted  as  the  limits  for  land- 
use  classes. 


Band  A 

Band  5 

Band  6 

Band  7 

Water 

0-127 

0-127 

0-127 

0-11 

Urban 

33-127 

17-127 

0-127 

8-63 

Woods 

0-127 

0-1A 

0-127 

8-63 

Open  Space 

0-32 

15-127 

0-127 

8-63 

Note:  Port  Allen  Vicinity  Classification,  Landsat 

1070-16073,  October  1972. 

These  class  ranges  were  used  to  test  each  pixel  in  the  Port  Allen  site 
in  an  automated  procedure.  When  a  pixel  was  within  the  class  limits  in 
all  four  bands,  it  was  assigned  that  class  code.  When  all  pixels  had 
been  tested  and  recorded  on  tape,  the  tape  was  used  to  drive  a  film 
writer  to  produce  a  unique  shade  of  gray  for  each  class,  making  a  picture 
of  the  land-use  classes  in  the  site.  Figure  11  is  such  a  picture  of 
Port  Allen  for  water,  forest,  and  other  classes.  Figure  11a  and  lib 
were  classified  with  100-m  grids  and  11c  and  lid  were  classified  with  a 
smaller  grid  (50  m) ,  which  is  nearer  the  original  pixel  size.  Figure  12 
separates  the  "other”  class  into  the  higher  reflectance  values  of  urban 
and  open  space  in  band  A. 

32.  This  classification  of  1972  data  was  used  in  the  enlarged 
site  to  include  Baton  Rouge  and  vicinity.  The  three-class  separation 
is  shown  in  Figure  13  and  the  four-class  separation  in  Figure  1A. 


Classification  Procedures  of  the  Enlarged  Site, 
Baton  Rouge  and  Vicinity,  1977  Data 


33.  Aerial  photography  over  the  Baton  Rouge  area,  flown  in  Febru¬ 
ary  197A,  was  located  at  EROS.  This  was  the  most  recent  photography 
found  for  use  in  terrain  sample  location  when  this  study  was  begun. 
Figure  15  shows  the  aerial  photograph  of  the  enlarged  Baton  Rouge  site 
with  the  accordance  transfer  points  identified:  A  is  the  midpoint  of  a 
bridge  over  a  portion  of  a  lake;  B  is  the  midpoint  of  a  tributary  enter¬ 
ing  the  Mississippi  River;  and  C  is  the  midpoint  of  a  tributary  entering 
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a.  Photography  -  February  1972 
100-m  grid 


b.  Landsat  -  October  1972 


100-m  grid 

Figure  11.  Water  (black),  woods  (white), 
and  other  (light  gray)  land-use  classes  at 
the  Port  Allen  site  (Sheet  1  of  3) 


a.  Photography  -  February  1972 
100-m  grid 


b.  Landsat  -  October  1972 
100-m  grid 


Figure  12.  Water  (black),  woods  (white),  open 
space  (light  gray),  and  urban  (dark  gray) 
land-use  classes  at  the  Port  Allen  site 
(Sheet  1  of  3) 


Landsat  -  October  1972 


igure 


50-m  grid 

12.  (Sheet  3  of  3) 


-c 


1 


a.  Photography  -  February  1972 


Figure  14.  Water  (black),  woods  (white),  open  space  (light  gray),  and 
urban  (dark  gray)  land-use  classes  at  the  Baton  Rouge  site  (Continued) 


the  Intracoastal  Waterway.  Figure  16  shows  the  locations  of  land-use 
sample  sites. 

Land-use  sample  locations  and  data 


34.  Twenty-seven  land-use  sample  sites  were  selected  from  the 
February  1974  aerial  photograph  (Figure  16).  Four  samples  (1-4)  were 


Figure  16.  Location  of  27  land-use  sample  sites  (1-27)  at  the 
Baton  Rouge  site,  February  1974  photography 


selected  in  what  appeared  to  be  heavily  wooded  areas;  ten  (5-14)  were  in 
the  city  of  Baton  Rouge;  seven  (15-21)  were  in  open  spaces,  mostly  sugar¬ 
cane  fields;  and  six  were  in  water  (22-27).  Equations  1  and  2  were  used 
with  the  appropriate  constants  from  Table  2  and  distances  of  the  sample 
locations  from  A  (Figure  15)  to  determine  the  corresponding  location 
(pixel  and  scan)  in  the  Landsat  data. 

35.  Reflectance  values  of  the  terrain  samples  and  their  eight  sur¬ 
rounding  terrain  pixels  were  located  in  a  digital  listing  of  the  Landsat 
data.  The  3  by  3  arrays  of  data  from  sample  1  (woods)  are  listed  below. 

Band  4  Band  5  Band  6  Band  7 

12  12  14  13  14  14  19  19  21  10  10  11 

13  11  13  13  15  13  18  21  21  10  10  11 

13  13  13  14  14  14  18  20  18  10  10  10 

Spectral  data  analysis 

36.  Water  and  woods.  Since  the  first  objective  of  the  study  was 
to  test  the  classification  procedure  to  identify  water  and  woods  at  the 
Baton  Rouge  site,  only  those  terrain  types  were  considered  at  this  point. 
The  mean  and  standard  deviation  were  computed  for  each  3  by  3  data  array 
to  determine  if  the  amples  were  representative  of  the  terrain  type  from 
which  they  were  selected  on  the  aerial  photo.  The  four  samples  from 
woods  (1-4  in  Table  6)  had  standard  deviations  less  than  1.5.  The  water 
samples  (22-27  in  Table  6)  contained  two  with  deviations  much  higher, 
samples  24  and  26.  Sample  24  was  from  a  lake  in  Baton  Rouge,  and  sam¬ 
ple  26  was  from  the  Intracoastal  Waterway.  Both  of  these  water  bodies 
are  so  narrow  that  the  pixel  values  are  mixtures  of  reflectance  from 
water  and  nonwater.  Having  selected  the  samples  in  the  widest  parts  of 
the  water  bodies,  there  seemed  to  be  only  two  choices:  delete  the  sam¬ 
ples  from  consideration,  or  reduce  the  array  by  eliminating  those  pixels 
suspected  of  having  reflectance  mixed  with  that  of  surrounding  terrain 
other  than  water.  The  samples  were  kept  and  band  7,  which  absorbs  almost 
all  reflectance  from  water,  was  analyzed.  Those  pixels  with  values 
greatly  different  from  the  other  band  7  values  were  eliminated.  The  same 
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pixels  were 

then  eliminated 

from  the 

other  arrays 

as  shown 

by 

the 

underlined 

elements 

in  the  following 

arrays : 

Sample  24: 

Band 

4 

Band 

5 

Band 

6 

Band 

7 

12  14 

14 

12  12 

13 

8  8 

15 

1 

1 

5 

12  13 

13 

13  12 

12 

8  8 

9 

0 

0 

2 

13  13 

13 

12  12 

11 

8  8 

8 

0 

0 

] 

Sample  26: 

Band 

4 

Band 

5 

Band 

6 

Band 

7 

17  15 

13 

24  20 

15 

18  18 

20 

4 

5 

9 

17  17 

13 

23  21 

15 

20  16 

21 

5 

3 

10 

14  17 

14 

17  22 

16 

22  16 

24 

9 

4 

10 

37.  The  means  and  standard  deviations  were  recomputed  for  the  two 
samples  and  listed  in  Table  6  as  revisions.  All  water  samples  now  met 
the  1.5  criterion  for  representativeness  except  sample  26,  which  had  a 
standard  deviation  as  high  as  1.67.  A  screening  window  of  1.75  from  the 
mean  was  accepted,  and  the  range  was  computed  for  each  band  of  all  re¬ 
tained  samples  (Table  7).  The  values  were  rounded  to  the  nearest  integer 
to  be  compatible  with  the  CCT's  values. 

38.  Next,  the  ten  woods  and  water  samples  were  tested  by  the  clus¬ 
tering  program  for  inclusiveness  of  terrain  type  and  exclusiveness  between 
types.  A  diagonalization  process*  grouped  those  samples  exhibiting  simi¬ 
larity  (i.e.,  samples  with  reflectance  from  all  bands  falling  within  the 
same  variance  windows).  A  similarity  matrix  plot  was  prepared  by 
placing  a  black  dot  where  the  vertical  and  horizontal  elements  agreed. 

The  matrix  in  Figure  17a  showed  that  a  variance  of  +1.75  would  not  form 
clusters  either  of  all  woods  or  of  all  water  samples.  The  variance 
window  was  increased  by  increments  of  0.25  until  2.25  included  all  woods 
samples  in  a  cluster  and  all  water  samples  except  sample  24  in  another 


cluster  (Figure  17b).  The  variance  window  was  increased  as  before  until 
3.75  was  reached  (Figure  17c).  Sample  24  was  then  similar  to  only  one 
other  water  sample,  and  the  two  clusters  were  no  longer  exclusive.  Sam¬ 
ple  24  was  shifted  outside  the  water  samples  cluster  and  a  diagonaliza- 
tion  at  variance  2.25  was  produced  as  shown  in  Figure  17d.  Since  sam¬ 
ple  24  did  not  encroacii  upon  the  characterization  of  woods,  it  was  kept 
in  the  analysis  as  a  water  sample. 
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a.  Similarity  matrix 
variance  =  1.75 


b.  Similarity  matrix 
variance  =  2.25 
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c.  Similarity  matrix 
variance  =  3.75 


d.  Diagonalization  by  cluster 
variance  =  2.25 


Figure  17.  Clusterinj  and  diagonalizing  of  wood  and  water  samples, 

Landsat  1977  data 


39.  The  variance  windows  of  each  sample  were  revised  to  +2.25  and 
new  class  variances  were  determined  (Table  7).  It  can  be  seen  that  the 
two  classes  can  be  separated  from  each  other  in  band  7;  in  fact,  none  of 
the  samples  except  water  in  Table  6  have  reflectance  values  in  the  0  to 
7  range.  Reflectance  values  in  band  7  from  0  to  7  characterize  water 
and  from  8  to  63  characterize  nonwater.  It  appears  that  woods  can  be 
discriminated  from  other  nonwater  pixels  by  band  5  with  a  class  window 
of  12  to  16.  A  check  of  all  sites  in  Table  6  shows  no  site  other  than 
woods  or  water  with  values  below  16  for  band  5.  Urban  sample  10  has  a 
mean  of  15,  but  it  is  a  sample  from  woods  within  the  urban  area.  By 
stretching  all  class  ranges  to  be  inclusive  in  order  to  classify  all 
pixels,  the  following  classification  was  developed  for  the  Baton  Rouge 


and  a  land-use 

map  produced 

(Figure  18) . 

A  visual 

comparison  of 

Band  4 

Band  5 

Band  6 

Band  7 

Water 

0-127 

0-127 

0-127 

0-7 

Woods 

0-127 

0-16 

0-127 

8-63 

Nonwater 

Nonwoods 

0-127 

17-127 

0-127 

8-63 

Figure  18  with  Figure  16  indicates  that  a  good  separation  of  water  and 
woods  from  each  other  and  from  other  land  uses  was  accomplished.  A  more 
complete  comparison  will  be  made  later  in  Part  III. 

40.  Urban  and  open  areas.  Another  objective  of  this  study  was  to 
explore  the  possibility  of  adding  urban  and  open  areas  to  the  classifi¬ 
cation  procedure.  Delineating  urban  areas  is  difficult  because  the 
reflectance  from  such  areas  is  not  unique  but  resembles  that  of  other 
land-use  types  in  a  varied  array  over  the  urban  area.  Samples  5-21  in 
Table  6  are  from  a  variety  of  land-use  types  in  Baton  Rouge  and  from 
open  fields.  It  can  be  seen  immediately  that  a  screening  window  of  1.5 
will  not  work  and  increasing  the  window  to  3.0  will  eliminate  almost 
half  of  samples  5-21.  A  review  of  the  reflectance  data  from  Baton  Rouge 
indicated  that  it  would  be  almost  impossible  to  get  homogeneous  3-  by  3- 
reflectance  samples  from  urban  areas.  Table  8  shows  7  by  7  arrays  of 
pixels  from  an  area  in  Baton  Rouge.  The  values  are  almost  a  checkerboard 


Figure  18.  Classification  of  water  (black),  woods  (gray),  and  other 
(white)  land-use  classes  at  the  Baton  Rouge  site,  1977 
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of  high  and  low  values.  This  condition  appeared  to  be  consistent  over 
the  whole  city  area.  High  and  low  variance  was  also  observed  in  the 
data  from  fields  where  bare,  vegetated,  wet,  and  dry  areas  were 
intermixed . 

41.  Ignoring  the  screening  window,  a  clustering  procedure  was 
begun  with  all  samples  except  water.  Water  was  not  included  because  it 
can  be  discriminated  from  all  other  samples  by  band  7.  With  a  variance 
of  1.75  (Figure  19a),  a  good  separation  could  be  seen  of  woods  (sam¬ 
ples  1-4)  from  open  (16-21)  samples.  There  seemed  to  be  some  grouping 
of  the  various  urban  samples.  Sample  5  showed  similarity  to  both  open 
space  and  to  urban  samples;  therefore,  it  was  eliminated.  The  variance 
was  extended  to  2.00,  and  the  clusters  began  to  lose  their  exclusiveness 
(Figure  19b).  Sample  10  showed  a  similarity  with  woods  and  not  urban. 
This  could  be  expected  since  it  was  a  sample  from  woods  in  the  urban 
area.  For  fear  of  distorting  the  classification  for  woods,  it  was  elimi¬ 
nated  from  the  procedure.  Samples  9  and  15  showed  similarity  to  each 
other  but  to  neither  of  the  types  from  which  they  were  drawn,  urban  and 
open,  respectively.  They  were  both  eliminated,  and  the  surviving  sam¬ 
ples  are  shown  in  a  diagonalized  cluster  array  with  a  variance  of  1.75 

in  Figure  20. 

42.  The  surviving  samples  with  variance  windows  of  1.75  are  listed 
in  Table  9  with  the  class  ranges.  These  ranges  were  plotted  (Figure  21) 
to  determine  unique  ranges  for  urban  and  open  space  to  be  added  to  the 
successful  class  ranges  of  woods  and  water.  Again  there  seemed  to  be  no 
way  to  separate  open  from  urban  except  for  those  urban  samples  with  very 
high  reflectance  such  as  samples  6,  7,  13,  and  14  (industrial  areas, 
large  flat  roofs,  parking  lots,  and  wide  concrete  thoroughfares).  Since 
these  terrain  types  are  found  in  urban  areas  and  not  in  open  spaces  such 
as  fields,  they  could  be  the  "locators"  for  the  urban  areas  and  perhaps 
even  define,  to  some  degree,  the  extent  or  borders  of  the  urban  areas. 

In  the  spaces  between  the  high  reflectance  pixels  there  are  many  pixels 
with  varying  lower  reflectance  that  complete  the  urban  area.  The  class 
range  for  urban  was  then  redefined  as  high  reflectance  in  band  4  (shown 
by  *  in  Table  9  and  indicated  by  a  double  bar  in  Figure  21).  Band  4 
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Figure  21.  February  1977  Landsat  range  of  reflectance  by 
terrain  class,  Baton  Rouge  site.  (Double  bar  defines  the 
accepted  urban  range) 
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shows  the  widest  separation  of  open  and  the  new  urban  range,  so  the 
value  halfway  between  them  rounded  to  the  nearest  integer  (23  or  24) 
could  be  accepted  as  the  lower  urban  class  limit.  After  several  runs, 
23  was  kept  as  the  lower  urban  limit. 

43.  By  increasing  the  limits  to  include  all  pixels,  the  classifi 
cation  becomes: 


Band  4 

Band  5 

Band  6 

Band 

Water 

0-127 

0-127 

0-127 

0-7 

Urban 

23-127 

17-127 

0-127 

8-63 

Woods 

0-127 

0-16 

0-127 

8-63 

Open 

0-22 

17-127 

0-127 

8-63 

Each  pixel  in  the  data  array  was  tested  for  compliance  with  land-use 
windows  in  order  (first,  water;  and  then  urban,  woods,  and  open  space) 
until  it  qualified  as  one  of  the  classes.  When  a  pixel  qualified  for  a 
land-use  class,  no  more  testing  was  performed  on  it  and  the  next  pixel 
in  the  scan  was  tested  and  so  on  to  the  next  scan  until  all  had  been 
tested.  All  pixels  will  fall  within  some  class  since  the  classification 
is  inclusive. 

44.  Figure  22  shows  the  results  of  this  classification.  The  green 
areas  indicate  woods;  blue,  water;  yellow,  open  spaces;  and  red  indi¬ 
cates  high  reflectance  surfaces,  most  of  which  are  located  within  urban 
areas.  The  high  reflectance  areas  outside  of  urban  areas  are  mostly 
from  highways  (especially  intersections),  oil  fields  and  refineries,  and 
clusters  of  buildings.  Some  red  areas  appear  to  be  from  sand  (in  the 
bend  of  the  river  in  the  top  center  of  the  image)  and  from  highly  re¬ 
flecting  areas  in  the  open  spaces,  e.g.,  bare  dry  soil.  These  must  be 
accepted  as  misclassifications  of  urban  areas  at  present. 

45.  Figure  23  is  a  picture  of  the  entire  1977  Landsat  scene, 
2768-15410,  using  the  same  classification  as  for  the  Baton  Rouge  area 
with  the  same  color  identifications.  Part  of  Lake  Pontchartrain  is  in 
the  center  of  the  east  edge.  The  Mississippi  River  runs  from  top  left  to 
just  below  Pontchartrain  on  the  right.  The  other  large  river,  westward 
of  the  Mississippi,  is  the  Atchafalaya.  Water  and  woods  are  well 
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Figure  23.  Classification  of  full  scene  Landsat  2768-15410  using 
classification  for  the  Baton  Rouge  site.  Green  areas  indicate 
woods;  blue,  water;  yellow,  open  space;  and  red,  urban 
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defined.  The  large  Atchafalaya  Basin  shows  almost  entirely  as  water, 
woods,  and  a  mixture  of  water  and  woods  with  very  little  open  space. 

Open  spaces  are  concentrated  around  the  Mississippi  River,  other  rivers, 
and  highways.  Clusters  and  strings  of  red  fall  where  cities,  commu¬ 
nities,  river  banks  (possibly  due  to  sand),  and  other  highly  reflecting 
terrain  are  located. 

Classification  Without  the  Guidance  of  Land-Use  Samples 


46.  Suppose  that  aerial  photography  or  other  terrain  information 
were  not  available  for  guided  classification.  Some  other  technique 
would  be  needed  for  reflectance  separation  of  terrain  classes.  A  pre¬ 
liminary  investigation  was  undertaken  with  the  same  Baton  Rouge  area 
Landsat  data  (28  February  1977)  to  find  those  pixel  signatures  that 
could  be  related  to  terrain  classes.  The  classification  procedure  was 
modified  to  substitute  a  computer  analysis  for  the  aerial  photography 
interpretation  of  terrain  classes  to  identify  a  reflectance  range  for 
each  band. 

47.  The  computer  analysis  was  made  in  the  following  steps: 

a.  Arrays  of  3  by  3  pixels  were  examined  by  the  computer  for 
similarity  with  a  given  standard  deviation  of  1.5.  The 
standard  deviation  is  a  variable  input  to  the  automated 
procedure . 

b.  The  means  of  those  arrays  with  similarity  were  listed  by 
homogeneity  within  each  band  beginning  with  band  7  and 
proceeding  to  band  4. 

c.  A  frequency  list  was  compiled  of  signatures  composed  of 
array  means  . 

d.  A  cluster  analysis  was  performed  on  those  signatures 
occurring  five  or  more  times. 

e.  The  diagona lization  procedure  was  applied  to  the  developed 
clusters . 

f.  Class  ranges  were  established  for  unique  clusters. 

g.  A  classification  picture  was  written  and  printed  for  the 
Baton  Rouge  site. 

48.  The  Baton  Rouge  site  contained  14,338  arrays  (320  scans  of 
403  pixels  each)  of  3  by  3  pixels  each.  Of  these,  4062  showed 
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similarity  with  a  1.5  standard  deviation.  Higher  standard  deviations 
could  have  been  selected,  which  would  have  been  more  inclusive,  or  lower 
for  more  restriction.  An  example  of  the  homogeneity  listing  is  presented 
in  Table  10,  and  an  example  of  the  frequency  listing  is  presented  in 
Table  11.  These  are  partial  listings.  As  can  be  seen  in  Table  10,  first 
the  signatures  with  band  7  mean  reflectances  from  lowest  to  highest  were 
listed,  with  sorting  continued  in  each  band  of  the  signature.  Band  6 
values  were  sorted  from  lowest  to  highest  in  each  unique  reflectance 
value  of  band  7,  and  so  on  until  all  signatures  with  homogeneity  had  been 
listed.  In  Table  11,  the  highest  number  of  any  one  signature  is  290  for 
signatures  13,  15,  21,  and  11  (reflectance  for  bands  4,  5,  6,  and  7, 
respectively) . 

49.  One  hundred  twenty-two  of  the  signatures  occurred  five  times 
or  more.  The  standard  deviation  was  reduced  to  1.3  and  84  signatures 


were  found 

occurring  five  or 

more  times 

Since 

1 . 3  may  not 

have  been 

out  of  the 

noise  level,  the 

signatures 

from  the 

1.5  standard 

deviation 

analysis  were 

used  in  the  clustering  procedure. 

50. 

To 

qualify  as  a  member  of  a 

cluster , 

a  signature 

must  exhibit 

similarity 

in 

all  bands  to  every  other 

member  in 

the  cluster 

Two  of  the 

one  hundred 

twenty-two  could 

not  qualify  in  this 

requirement 

and  were 

eliminated . 

Ten  separations 

were  found 

with  the 

clustering 

procedures  . 

The  ranges , 

lowest  to  the  highest  for  each  band 

in  a  cluster 

,  are  shown 

below: 

Number  of 

Range  of  Reflectance  Values 

Cluster 

Signatures 

Band  4 

Band  5 

Band  6 

Band  7 

1 

26 

17-19 

20-22 

12-14 

1-2 

2 

36 

12-14 

14-16 

19-22 

10-12 

3 

10 

12-14 

13-16 

19-21 

9-11 

4 

18 

13-15 

14-16 

23-25 

11-12 

5 

14 

19-21 

25-27 

28-30 

11-13 

6 

6 

19-20 

23-25 

27-29 

11-12 

7 

4 

14-15 

17 

23-25 

11-13 

6 

4 

12-13 

13-14 

17-19 

8-10 

9 

1 

21 

27 

31 

13 

10 

1 

17 

21 

25 

11 
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51.  From  this  point,  classification  was  no  longer  unguided.  Prior 
knowledge  of  reflectance  values  associated  with  land-use  classes  and  the 
grayscale  pictures  of  each  band  (Figure  6)  were  used  to  identify  clus¬ 
ters  by  land-use  classes.  Low  values  in  Sand  7  had  been  used  to  sepa¬ 
rate  water  from  other  reflectance  values;  therefore,  cluster  1  was  iden¬ 
tified  as  water.  Low  values  in  band  5  had  been  used  to  successfully  sep¬ 
arate  woods  from  other  reflectance  in  that  band;*  therefore,  clusters  2, 
3,  A,  and  8  were  combined  for  woods  identification: 


Band  A 

Band  5 

Band  6 

Band  7 

Cluster 

Water 

Woods 

17-19 

12-15 

20-22 

13-16 

12-1A 

17-25 

1-2 

8-12 

1 

2, 3, A, 8 

52.  A  look  at  Figure  6  indicates  that  band  A  shows  the  greatest 
contrast  between  pixels  in  Baton  Rouge  and  those  in  open  areas.  A  histo¬ 
gram  of  reflectance  values  indicated  that  the  similarity  signatures  did 
not  include  the  highest  reflectance  values.  Apparently  the  3  by  3  matri¬ 
ces  containing  the  high  values  of  band  A  did  not  meet  the  homogeneity 
criterion.  Upon  examination  of  the  data,  this  proved  to  be  true.  This 
was  understandable  since  it  is  probable  that  the  higher  values  came  from 
highly  reflective  surfaces  such  as  buildings,  hard-surfaced  highways, 
large  flat  rooftops,  and  parking  areas.  These  surfaces  are  complex  in 
arrangement  and  varying  in  slope  and  proximity,  i.e.,  not  in  a  homoge¬ 
neous  arrangement  such  as  water  or  woods.  The  remaining  similarity  clus¬ 
ters,  with  the  lower  values  in  band  A,  should  be  representative  of  open 
spaces  such  as  space  occupied  by  sugarcane  fields,  pastures,  and  crop¬ 
lands;  in  fact,  it  was  thought  they  must  be  associated  with  everything 


not 

water,  wood,  or  urban. 

Combining 

clusters  5,  6, 

7 ,  9 ,  and  10 

into 

one 

range  and 

assigning  the 

higher  reflectance  values 

in  band  A  to 

urban, 

the 

following 

associations 

were  made: 

Band  A 

Band  5 

Band  6 

Band  7 

Open 

1A-21 

17-27 

23-31 

11-13 

Urban 

22-127 

0-127 

0-127 

0-63 

*  Ibid. 
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53.  Figure  24  shows  plots  of  the  cluster  groupings  and  reflec¬ 
tance  values  from  paragraph  52.  Another  consideration  was  given  to 
those  higher  reflectance  values  not  appearing  in  the  similarity  signa¬ 
tures.  It  was  supposed  that  alder  residential  areas  in  Baton  Rouge 
would  have  trees,  which  would  cause  pixels  from  those  areas  to  be 
classed  as  woods  rather  than  urban,  but  would  have  higher  reflectance  in 
bands  4  and  6.  The  following  separation  was  made:  if  pixels  classed  as 
woods  by  band  5  had  higher  band  6  reflectance  than  those  included  in  the 
list  of  homogeneous  signatures,  then  they  were  classed  as  other  urban. 

54.  When  the  range  limits  were  increased  to  be  inclusive  of  all 
pixels,  the  following  classification  was  made: 


Band  4 

Band  5 

Band  6 

Band  7 

Water 

0-127 

0-127 

0-127 

0-7 

Woods 

0-127 

0-16 

0-25 

8-63 

Other  (urban) 

0-127 

0-16 

26-127 

8-63 

Open 

0-21 

17-127 

0-127 

8  63 

Urban 

22-127 

17-127 

0-127 

8-63 

Figure  25  shows  the  resulting  classification.  While  a  close  look  at  Fig 
ure  25  shows  some  of  the  other  urban  class  in  the  urban  area,  more  of  it 
appears  in  what  should  have  been  open  space.  No  further  study  was  made 
in  this  direction;  however,  a  comparison  of  Figure  25  with  Figure  15, 
the  guided  classification,  indicates  that  water  and  woods  can  be  classi¬ 
fied  very  well  without  guidance  from  terrain  samples. 
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Figure  24.  Reflectance  ranges  grouped  from  unguided  cluster  analysis 


Figure  25.  Baton  Rouge  area  classified  without  sample  guidance,  2768- 
15410.  Green  areas  indicate  woods;  blue,  water;  yellow,  open  space; 
brown,  other  urban;  and  red,  urban 


PART  III:  EVALUATION 


Evaluation  Procedure 

55.  Evaluations  of  land-use  classification  from  Landsat  data  were 
made  to  test  the  procedure  (in  areas  other  than  the  Satartia  test  site 
used  by  Struve,  Grabau,  and  West  (1977b)*)  for  classifying  water, 
forested,  and  other  land-use  classes.  The  evaluation  was  made  by  com¬ 
paring  the  Landsat  classified  maps  to  photo- interpreted  land-use  maps  of 
the  same  site  area.  Both  maps  were  rectified  to  UTM  coordinates  on  the 
same  scale  topographic  maps  and  gridded  to  the  same  cell  size  for  a  com¬ 
parison  of  areas  in  each  land-use  class. 

Photo- interpreted 
land-use  maps  (UTM) 

56.  Land-use  classes  (water,  woods,  urban,  and  open  space)  were 
manually  interpreted  from  aerial  photography  by  outlining  the  areas  on 
transparent  overlays  inside  a  rectangular  border  (area  of  interest).  The 
corners  of  the  area  of  interest  were  identified  by  UTM  coordinates  in 
thousands  of  metres.  The  class  areas  (patches)  on  the  overlay  were  then 
digitized;  i.e.,  their  borders  were  defined  in  terms  of  digitizer  coordi¬ 
nates  that  were  referenced  to  the  UTM  coordinates  of  the  upper  left  cor¬ 
ner.  The  digitized  data  were  stored  on  magnetic  tape  and  later  fed  into 
a  computer  program  that  transposed  the  patch  descriptive  data  into  a  rec¬ 
tangular  grid  system  by  coding  each  space  in  the  grid  system  with  the 
land-use  code  of  the  patch  containing  the  space.  The  grid  size,  an  op¬ 
tion  in  the  computer  program,  was  chosen  as  100  by  100  m  for  comparison 
with  the  Landsat  classification  maps.  The  gridded  data  were  stored  on 
magnetic  tape  and  fed  into  the  film  writer  to  produce  land-use  maps  (Fig¬ 
ures  13a,  14a). 

Landsat  land-use  maps  (UTM) 

57.  A  Landsat  land-use  map  for  the  same  area  of  interest  (de¬ 
scribed  by  the  same  UTM  corner  coordinates)  was  produced  by  the  follow¬ 
ing  steps. 

*  Ibid. 
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a.  Classification  of  the  full  Landsat  scene. 

(1)  A  computer  program  was  developed  to  join  the  WES  Land- 
sat  tapes  (west  and  east  halves)  of  the  Landsat  scene. 

(2)  The  full  Landsat  scene  was  classified  using  class 
reflectance  ranges  as  discussed  in  Part  II,  e.g.. 

Port  Allen  site  (paragraph  31). 

b.  Rectification  of  the  site  land-use  map  to  UTM  coordinates. 

(1)  Five  clearly  defined  points  (e.g.  intersections  of 
major  highways)  in  the  vicinity  of  the  selected  site 
were  located  on  a  topographic  map  and  identified  by 
their  UTM  coordinates. 

(2)  The  corresponding  points  were  located  on  the  Landsat 
band  5  image,  and  50-  by  50-pixel  data  arrays  centered 
on  each  point  were  extracted  from  the  reflectance  data. 

(3)  Using  the  high  reflectance  values  in  band  5  to  iden¬ 
tify  the  transfer  points,  the  pixel  and  scan  of  each 
of  the  five  points  were  determined. 

(4)  The  two  sets  of  identifying  data,  UTM  coordinates  and 
corresponding  pixels  and  scans,  were  fed  into  a  com¬ 
puter  program  that  computed  the  angle  of  rotation, 
the  width  and  length  of  the  transposed  pixel,  and  the 
UTM  coordinates  of  the  first  pixel  and  scan  of  the 
transposed  map.  These  constants  were  used  to  tran¬ 
spose  any  array  of  extracted  data  in  the  locale  of 
these  transfer  points  to  UTM  coordinates.  This  resam¬ 
pling  technique  was  applied  to  unstretched  classified 
data,  i.e.,  the  third  and  twentieth  repeated  scans 
which  were  used  for  geometric  accordance  (para¬ 
graph  26)  were  removed.  The  resampling  procedure 
also  includes  the  option  to  change  the  pixel  size  in 
the  transposed  map  to  selected  width  and  length.  In 
this  study,  100-  by  100-m  pixels  were  selected  and 
the  transposed  map  data  stored  on  magnetic  tape. 

Accordance  of  Landsat  map 
with  p hoto -interpreted  map 

58.  The  final  maps  were  then  prepared  with  the  photo-interpreted 
map  as  the  base  map  or  standard  for  comparison.  The  following  steps 
were  used: 

a.  The  Landsat  classified  data,  transposed  to  UTM  coordi¬ 
nates,  were  overlaid  by  the  data  of  the  aerial  photograph 
land-use  map,  which  was  used  as  a  mask  to  bring  the  two 
maps  into  area  agreement. 

b.  The  resulting  Landsat  classified  data  were  fed  into  the 
film  writer  and  a  land-use,  gray-shaded  map  was  written. 
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A  site  in  the  Atchafalaya  Basin  (Figure  26)  was  classi¬ 
fied  and  rectified  by  this  procedure  (Figure  27).  Fig¬ 
ure  28  is  a  photo-interpreted  land-use  map  of  the  same 
area  of  interest.  The  two  maps  now  covered  the  same  area 
with  the  same  .scale  and  were  rectified  to  the  same  UTM 
coordinates . 

Comparison  of  Landsat  map  with 
photo-interpreted  map  by  total  area 

59.  A  computer  program  was  written  to  sum  the  number  of  pixels 
(grid  spaces)  for  each  class  in  each  map,  to  convert  this  sum  to  acreage 
or  other  area  measurement,  and  to  compute  the  percent  of  each  class  in 
the  total  site.  The  percent  of  the  Landsat  classes  in  the  total  site 
were  compared  to  those  derived  from  the  photo-interpreted  procedure, 
which  were  accepted  as  the  true  percents.  The  difference  in  percents 
expresses  the  error  in  classification  by  Landsat.  Table  12  shows  the 
results  of  this  procedure  for  the  Port  Allen  and  Baton  Rouge  sites  and 
for  Satartia,  the  site  used  in  the  development  of  the  classification 
procedure  described  by  Struve,  Grabau,  and  West  (1977b).* 

Automated  compari- 
s  o n  of  land-use  maps 

60.  The  comparison  of  percent  area  of  the  site  is  meaningful  but 
may  lead  to  false  conclusions  in  classification  accuracy.  The  total 
area  by  class  may  agree,  but  the  positions  of  like  classified  grid 
spaces  may  be  drastically  different.  A  new  program  was  used  to  compare 
not  only  the  total  number  of  correctly  identified  grid  spaces  (area), 
but  also  to  compare  the  two  maps  grid-by-grid  in  the  area  of  interest 
common  to  both  maps.  The  program  will  accept  maps  of  different  size, 
orientation,  and  grid-space  dimensions.  The  required  input  to  the  pro¬ 
gram  consists  of  the  two  data  arrays,  a  set  of  control  points  describing 
common  ground  surface  points  in  terms  of  each  map's  x-  and  y-distances, 
land-use  codes  on  each  map,  patch  numbers  on  the  standard  map,  and  de¬ 
sired  land-use  code  comparisons.  In  this  study,  the  aerial  photography 
maps  were  the  accepted  standards. 

61 .  The  output  of  the  automated  procedure  is  a  series  of  tables 


*  Ibid. 
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Figure  26.  Reduced  Chicot  Lake  quadrangle  sheet  containing  a  site 

in  the  Atchafalaya  Basin 


V 


Figure  28.  Portion  of  Atchafalaya  Basin  classified  by  manual 
interpretation  of  aerial  photography.  (Black  represents 
water;  white,  forest;  and  gray,  open) 
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listing  the  control  point  locations,  the  comparison  of  x-  and  y- 
coordinates  of  the  two  maps,  the  number  of  correctly  classified  grid 
spaces,  the  land-use  classes  incorrectly  identified,  and  the  percent  of 
correctly  and  incorrectly  identified  spaces  by  total  areas,  by  patch 
areas,  and  by  class. 

62.  Table  13  is  an  example  of  the  control  points  of  a  standard 
map  and  another  map  (Landsat)  and  the  ratio  of  the  x-  and  y-coordinates 
of  the  two  maps.  In  this  example,  the  Landsat  map  area  of  interest  (171 
by  141  grid  spaces)  was  a  part  of  a  larger  (311  by  331  grid  spaces) 
aerial  photography  map. 

63.  Table  14  lists  the  differences  between  the  projected  and 
actual  distances  of  the  Landsat  control  points.  Points  5  and  6  are  as 
much  as  four  or  five  grid  spaces  different  from  other  points.  The 
greatest  misclassifications  will  likely  occur  in  the  map  areas  of  these 
control  points. 

64.  Table  15  lists  the  number  of  grid  spaces  in  the  Landsat  map 
that  were  correctly  identified  by  land-use  class.  Also  shown  are  the 
classes  assigned  to  the  incorrectly  identified  grid  spaces,  e.g.,  of  the 
2307  water  grid  spaces  on  the  standard  map,  1751  were  classed  correctly 
by  class  and  position;  94  were  classed  as  woods;  404  as  open  space;  and 
58  as  urban.  In  addition  to  the  1751  correctly  identified  water  grid 
spaces,  Landsat  incorrectly  identified  116  woods,  48  open  space,  and  20 
urban  spaces  as  water. 

65.  Figures  29-31  illustrate  the  comparison  technique.  Patches 
including  only  one  land-use  class  are  manually  drawn  with  their  class 
codes  on  an  overlay  to  the  photo  (Figure  29a).  The  boundaries  of  the 
patches  are  identified  in  terms  of  map  and  UTM  coordinates  by  a  semi- 
automated  digitizing  process  and  automatically  stored  on  magnetic  tape. 
The  digitized  information  is  used  in  an  automated  gridding  procedure 
(Figure  29b)  to  produce  a  map  (Figure  29c).  The  transformation  of  the 
Landsat  map  for  geometric  accordance  with  the  standard  map  is  performed 
sequentially  for  each  grid  space  in  the  standard  map.  The  class  code  of 
the  corresponding  Landsat  grid  space  is  transposed  to  a  new  map  coinci¬ 
dent  in  position  with  the  standard  map.  The  following  steps  are 
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performed  in  the  transformation  (Figure  30): 

a.  The  control  points  are  checked  for  the  nearest  four  that 
lie  in  different  quadrants  relative  to  the  grid  space 
(Figure  30a). 

b.  The  centroid  of  the  four  control  points  is  located  in 
both  standard  and  Landsat  maps  (Figures  30b  and  c) . 

£.  The  distance  and  direction  of  the  grid  space  from  the 
centroid  in  the  standard  map  are  modified  by  a  weighted 
square  of  the  distance  and  the  difference  in  scale  and 
rotation  of  the  Landsat  map  (Figure  30b). 

d.  The  modified  distance  and  direction  are  used  to  locate  the 
corresponding  grid  space  in  the  Landsat  map  (Figure  30c). 

e.  The  class  code  of  the  Landsat  grid  space  is  placed  in  the 
grid  space  of  the  transform  map  that  is  coincident  with 
the  grid  space  of  the  standard  map  (Figure  30d). 

f.  The  above  procedure  is  repeated  for  all  grid  spaces  in 
the  area  of  interest  of  the  standard  map. 

66.  The  class  of  each  grid  space  in  the  standard  map  is  compared 
to  the  class  of  the  corresponding  grid  space  in  the  transform  map  (Fig¬ 
ure  31).  Those  spaces  in  either  map  not  coincident  (e.g.,  line  1  in 
Figure  31)  are  excluded  from  the  summing  activity.  Lines  2,  3,  and  5 
join  coincident  spaces  with  the  same  classes  in  both  maps  and  will  be 
included  in  the  sum  of  the  correctly  identified  classes  in  the  patch  and 
in  the  map.  Line  4  joins  coincident  spaces  with  different  classes  and 
will  be  summed  accordingly,  i.e.,  summed  as  incorrect  by  class  in  the 
patch  and  in  the  map. 

67.  Table  16  is  a  partial  output  listing  of  the  number  of  spaces 
in  the  standard  map  patches,  the  correct  class,  and  classes  on  the  Land¬ 
sat  map  for  the  corresponding  patch.  Note  that  some  patch  numbers  are 
missing  from  the  list.  These  were  one-space  patches  and  were  eliminated 
from  the  correlation  procedure  by  an  island-removal  program  for  removal 
of  spurious  data.  The  program  tests  the  eight  pixels  surrounding  the 
isolated  space  and  assigns  it  to  the  class  most  frequently  occurring 
among  the  eight.  A  test  case  (3  classes)  of  before  and  after  island  re¬ 
moval  showed  only  a  small  change  in  the  summary  results,  from  91.696  per¬ 
cent  to  92.108  percent  for  a  gain  of  0.512  percent  correctly  classified. 

68.  Table  17  is  the  output  summary  of  the  class  correlation 
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patch-by-patch  of  the  two  maps.  The  list  is  arranged  by  sequence  of  the 
class  codes  and  by  size  within  the  class;  e.g.,  water  (code  1)  is  listed 
first  with  patch  13  first  within  the  class  (1559  grid  spaces).  The  total 
area  (standard  map),  the  number  of  correct  and  incorrect  classifications, 
and  the  percent  of  patch  and  percent  of  map  are  shown  for  correct  and 
incorrect  classifications. 

Order  in  which 
evaluations  were  made 

69.  The  October  1975  Satartia  site  results*  were  used  as  a  refer¬ 
ence  for  the  success  of  classification  by  area  in  this  study.  These 
were  from  a  three-class  separation:  water,  woods,  and  the  remainder  of 
the  site  in  the  class  called  "other."  The  first  evaluations  were  made 
for  the  Port  Allen  and  Baton  Rouge  three-class  separations  compared  to 
the  Satartia  site  results  for  area  percentage.  Next,  an  evaluation  was 
made  of  the  three-class  separations  using  area  percentage  and  position 
as  criteria.  Since  this  kind  of  result  was  not  made  by  Struve,  Grabau, 
and  West  (1977b),*  only  the  comparison  to  the  photo-interpreted  maps  was 
available  as  a  standard.  Following  the  evaluations  of  the  three-class 
separations,  the  four-class  separation  evaluations  were  made  by  area 
percentage  and  position. 

70.  Other  evaluations  were  then  made  of  (a)  the  difference  in 
class  ranges  of  reflectance,  (b)  the  effect  of  shape  and  size  of  class 
patches,  (c)  the  orientation  of  patches  to  the  scanning  of  Landsat  sen¬ 
sors,  (d)  thi  limitations  of  photo  interpretation  as  a  standard,  (e)  the 
limitations  of  the  Landsat  classification  and  change  detection,  and  (f) 
the  use  of  Landsat  to  identify  urban  areas. 

Water  and  Woods  Classification  Comparisons 

71.  The  comparisons  by  area  of  the  Landsat  classifications  with 
photo  interpretation  classifications  for  three  classes  are  given  in 
Table  12  for  Satartia  and  Baton  Rouge  sites.  The  two  sites  are  similar 
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in  climate,  soil,  and  physiography.  They  both  lie  in  the  flat,  low- 
lying  floodplain  of  the  Mississippi  River  with  soil  formations  of  the 
Pleistocene  age.  The  summers  are  long  and  hot.  Precipitation  is  the 
heaviest  in  the  winter  and  spring  months,  lighter  in  the  summer  months, 
with  an  average  annual  precipitation  of  53  and  54  in.,  respectively,  for 
Satartia  and  Baton  Rouge.  Satartia  lies  approximately  2  deg  almost  due 
north  of  Baton  Rouge.  Cotton  and  soybeans  are  major  crops  in  the  vicin¬ 
ity  of  Satartia.  The  major  crop  in  the  Baton  Rouge  vicinity  is  sugar¬ 
cane.  The  major  difference  in  land  use  of  the  two  areas  is  the  large 
urban  area  in  the  Baton  Rouge  site.  With  the  similarity  of  the  two 
sites,  a  close  comparison  in  classification  was  expected  for  the  data 
collected  in  October  for  both  sites.  All  of  the  guided  classifications 
were  within  2.2  percent  of  the  Satartia  results.  The  classification  re¬ 
sults  of  the  1972,  50-m  grid  for  Port  Allen  and  the  100-m  grid  for  Baton 
Rouge  were  closer  to  the  photo  interpretation  than  Satartia  as  shown  in 
Tables  12  and  18. 

72.  Comparison  by  position  decreases  the  percentage  of  correct 
classification  by  6  to  7  percent  for  1972  and  by  9.2  percent  for  1977 
(Table  19).  Pictures  of  the  Landsat  classification  and  the  photo- 
interpreted  standard  maps  for  the  Port  Allen  and  Baton  Rouge  sites  are 
shown  in  Figures  11,  13,  and  32.  The  map  data  for  these  pictures  were 
used  in  the  comparison  program  giving  the  three  class  comparisons  of 
Table  19  and  the  following  results: 

Percent  of  Site  Correct 


Port  Allen  1972, 

100  m 

91.7 

Port  Allen  1972, 

50  m 

91.5 

Baton  Rouge  1972 

92.1 

Baton  Rouge  1977 

85.9 

73.  Some  thought  should  be  given  to  the  selection  of  a  standard  by 
which  the  Landsat  classification  is  measured.  Difficulty  exists  in  ob¬ 
taining  photography  taken  at  the  proper  time  and  of  reliable  quality; 
the  alternative  is  to  bear  the  expense  of  having  photographic  missions 
flown.  Even  with  the  best  photography  at  the  desired  time,  there  are 


a.  Photography  -  October  1974 


Figure  32.  Water  (black),  woods  (white),  and  other  (gray) 
land-use  classes  at  the  Baton  Rouge  site  (Continued) 


still  other  factors  affecting  its  quality  as  a  standard,  such  as  the 
manual  interpretation  of  land-use  classes,  the  production  of  the  class 
patch  maps  for  the  digitizing  function,  the  digitizing  itself  (whether 
with  a  hand-held  sensor  or  with  an  electronic  scanner),  the  conversion 
of  digitized  data  on  tape  to  a  machine-compatible  format,  and,  finally, 
the  interpretation  of  the  digitized  data  to  an  array  representing  the 
land  surface  by  a  gridding  procedure.  There  are  many  places  for  human 
error,  perhaps  more  than  in  the  classification  of  the  Landsat  data. 


Urban  and  Other  Classification  Comparisons 


74.  The  class  called  "other”  in  the  three-class  separation  in¬ 
cluded  urban  and  open  space.  In  evaluation  by  area  as  percent  of  the 
total  site,  the  percent  correct  remained  good  for  the  three-class  separa 
tion  when  using  the  new  automated  comparison  program.  Results  are  shown 
as  follows: 

Percent  of  Site  Correct 

Port  Allen  1972,  100  m  94.7 

Port  Allen  1972,  50  m  95.1 

Baton  Rouge  1972  84.7 

Baton  Rouge  1977  72.7 

75.  The  aerial  photography  used  for  the  1977  Landsat  classifica¬ 
tion  sample  identification  was  very  poor.  It  was  color  infrared, 
1:131,000  enlarged  to  1:65,500,  flown  in  February  1974,  but  it  was  the 
most  recent  found  at  the  beginning  of  the  study.  The  1972  photography 
was  of  much  better  quality. 

76.  One  of  the  biggest  factors  in  the  poor  results  of  the  four- 
class  separation  comparisons  is  the  interpreter's  criterion  that  indus¬ 
trial,  business,  and  all  residential  areas  (whether  new  or  old)  make  up 
the  urban  class.  The  Landsat  sensors  record  the  reflectance  as  it  is; 
an  older  residential  area  with  a  lot  of  trees  would  be  classed  as  woods 
even  though  it  is  next  to  an  industrial  area.  The  photo  interpreter 
calls  it  urban;  Landsat  classification  did  not  make  that  interpretation. 


The  second  biggest  factor  is  probably  the  similarity  of  integrated  re¬ 
flectance  from  a  mixture  of  entirely  different  land  use.  Many  reflec¬ 
tance  values  in  Baton  Rouge  from  residential  areas  are  from  tree  tops, 
rooftops,  and  paved  streets.  The  resultant  integration  is  reflectance 
higher  than  that  from  woods,  lower  than  that  from  highways  and  buildings, 
but  very  like  that  of  the  integration  of  reflectance  from  sugarcane 
fields,  bare  ground,  and  a  few  farm  buildings.  This  can  be  seen  in  the 
pictures  of  the  aerial  photography  and  Landsat  classifications  in  Fig¬ 
ures  12,  14,  and  33  and  by  referring  to  Figures  2  and  7. 

77.  The  comparison  of  the  four-class  separation  by  position  shows 
a  great  degradation  of  the  percent  of  the  site  correct.  The  results  are 
shown  below: 

Percent  of  Site  Correct 


Port  Allen  1972, 

100  m 

74.3 

Port  Allen  1972, 

50  m 

73.6 

Baton  Rouge  1972 

71.6 

Baton  Rouge  1977 

59.5 

78.  Another  factor  in  the  poor  correlation  of  the  four-class 
separation  using  the  new  comparison  by  position  may  be  the  difference  in 
the  gridding  procedures  for  photo  interpretation  and  Landsat  classifica¬ 
tion.  The  Landsat  array  was  rectified  to  UTM  coordinates  from  an  exist¬ 
ing  array  by  locating  common  points  in  the  two  data  sets  and  computing 
the  rectification  constants.  The  photo-interpreted  array  was  constructed 
from  digitized  patch  data  with  a  set  of  rules  governing  the  border  pixel 
selection.  It  was  brought  into  accordance  with  UTM  coordinates  manually 
from  a  topographic  map. 


Other  Evaluations 

Size  and  shape  of  pixel  and  patches 

79.  Reducing  the  size  of  the  grid  (paragraph  65)  seems  to  have 
very  little  effect  on  the  accuracy.  The  three  classes  by  area  comparison 
show  a  slight  improvement  as  shown  in  Tables  12  and  18.  However,  the 
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a.  Photography  -  October  1974 

Figure  33.  Water  (black),  woods  (white),  open  space  (light  gray), 
and  urban  (dark  gray)  land-use  classes  at  the  Baton  Rouge  site 

(Continued) 


area  and  position  comparison  shows  a  loss  of  accuracy. 

80.  The  size  and  shape  of  patches  appear  to  have  an  influence  on 
the  accuracy  of  classification.  In  Figure  13,  the  opportunities  for  mis¬ 
calculations  occur  along  the  interface  of  the  "other"  class  with  water 
and  woods.  The  number  of  pixels  along  the  interface  is  small  compared 

to  the  total  "other"  class,  which  means  a  small  percent  of  misclassif ica- 
tions.  Patches  with  many  irregularities  along  the  border  resulting  in  a 
large  number  of  pixels  along  the  interface  compared  to  the  number  within 
the  border  of  the  patch  have  more  opportunity  for  misclassif ication. 

This  may  account  for  the  lower  accuracy  of  the  50-m  pixel  classification 
of  Port  Allen  for  the  four-class  separation.  Compare  the  photo  interpre¬ 
tations  in  Figures  11  and  12. 

81.  The  effect  of  shape  can  be  seen  in  Figure  12.  In  Figure  12a, 
the  photo  interpretation  shows  three  one-pixel-wide  roads  (darkest 
shade):  one  in  the  north-south  direction  on  the  upper  west  side  of  the 
scene,  another  in  the  east-west  direction  in  the  southern  part  of  the 
scene,  and  the  third  nearby  in  a  slightly  southeast  to  northwest  direc¬ 
tion.  The  Landsat  scanner  will  cross  this  scene  at  9.76  deg  east  of 
north  (Table  2).  Figure  12b  shows  that  the  long  narrow  shape  of  the 
first  road  with  orientation  across  the  scan  has  almost  disappeared  as 
has  the  second  road.  A  good  portion  of  the  third  road  has  been  identi¬ 
fied.  This  road  is  closer  to  paralleling  the  scan  line  and,  even  though 
it  is  only  one  pixel  wide,  it  was  picked  up  and  classified  as  the  road 
patch.  In  Figure  12c,  the  north-south  road  is  two  pixels  wide  and  the 
Landsat  interpretation  retains  the  road  in  Figure  12d. 

Detection  of  changes 

82.  An  evaluation  was  made  of  the  ability  of  Landsat  to  detect 
change.  Two  good  examples  are  shown  between  the  1972  and  1977  Baton 
Rouge  sites.  In  the  upper  right  corner  of  the  site  in  Figure  2,  the  new 
interstate  highway  leading  out  of  Baton  Rouge  to  the  northeast  can  be 
seen  to  come  to  an  abrupt  end  at  a  wooded  patch  where  the  construction 
ended.  The  Landsat  classification  shows  the  wooded  patch  although  only 
part  of  the  highway  is  shown  (Figure  14).  The  1974  aerial  photography 
(Figure  15)  shows  the  road  under  construction  moving  into  the  wooded 
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patch,  and  the  1977  Landsat  shows  the  road  completed  through  the  westr 
side  of  the  wooded  patch  in  Figure  33  and  more  clearly  in  the  color  pic¬ 
tures,  Figures  22  and  23. 

83.  Another  example  of  change  detection  is  in  the  Port  Allen  site. 
In  the  1972  photographs,  north  of  Highway  10  and  west  of  Highway  1,  new 
subdivisions  with  few  trees  and  bare  ground  give  a  high  reflectance,  as 
can  be  seen  in  Figures  2,  7,  and  12.  It  is  speculated  that  by  1977 
these  subdivisions  had  older  trees  and  lawns,  causing  the  reflectance 
values  to  be  lower  and  to  be  classed  as  open  space.  In  Figure  33,  this 
area  retains  the  high  reflectance  of  major  streets,  but  has  lost  the 
high  reflectance  it  had  when  the  subdivision  was  new.  Again  Figures  22 
and  23  show  the  change  from  Figures  2  and  12. 

Limited  number  of  classes 

84.  A  look  at  the  class  ranges  of  each  site  revealed  that  both 
February  and  October  Satartia  scenes  had  class  ranges  for  woods  and  other 
classes  slightly  higher  than  the  Baton  Rouge  scenes.  Table  20  lists  all 
class  ranges  used  in  this  study.  The  decision  to  limit  the  number  of 
classes  to  three  and  four  and  to  force  the  inclusiveness  of  all  pixels 

in  the  site  results  in  a  loss  of  information  about  where  misclassifica- 
tion  is  most  likely  to  occur.  Perhaps  a  better  understanding  could  be 
obtained  (concerning  the  mixture  of  reflectance  values  as  one  class 
gradually  becomes  dominant  in  a  land  space)  by  leaving  the  ranges  un¬ 
stretched  and  studying  the  signatures  of  these  unclassified  pixels, 
rather  than  making  all  pixels  inclusive.  Or,  perhaps  relationships  be¬ 
tween  bands  in  the  mixed  reflectances  could  be  established. 

Photo  interpretation  limitations 

85.  Some  limitations  of  photo  interpretation  were  discovered. 

The  quality  of  available  aerial  photography  is  vital  to  the  reliability 
of  the  standard  map.  Manual  input  required  at  several  stages  of  develop¬ 
ing  a  photo-interpreted  standard  can  be  the  source  of  errors.  Two  of 
these  manual  inputs  are  the  land-use  interpretation  and  the  digitizing 
effort  to  place  the  data  on  magnetic  tapes  for  use  in  the  computer. 

86.  A  cursory  evaluation  was  made  of  the  classification  of  a  site 
in  the  Atchafalaya  Basin  that  was  on  the  Landsat  scenes  used  in  this 
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study.  The  area  is  shown  in  Figure  26,  which  is  a  reduction  of  the 
quadrangle  containing  the  site.  The  site  includes  large  swampy  areas 
containing  woods  and  water.  A  photo  interpretation  was  made  as  shown  in 
Figure  28.  The  interpreter  was  unable  to  show  all  the  small  patches  of 
water  and  approximated  total  areas  in  several  medium  size  patches.  Land- 
sat  classification,  shown  in  Figure  27,  was  able  to  distinguish  water 
not  visible  to  the  interpreter  in  the  photography.  This  could  have  been 
caused  either  by  accumulation  of  water  in  the  low-lying,  flat  areas  be¬ 
tween  the  time  of  the  photograph  and  the  Landsat  overpass,  or  by  the 
lack  of  experience  of  the  person  doing  the  photo  interpretation.  The 
Landsat  classification  gave  a  more  accurate  picture  of  the  area  in 
digital  data  that  can  be  summed  for  more  accurate  acreage  determination. 
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PART  IV:  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 


87.  The  classification  procedure  described  by  Struve,  Grabau,  and 
West  (1977b)*  and  developed  with  data  from  the  Satartia  site  works 
equally  as  well  for  the  Baton  Rouge  site  and  the  Satartia  site  to  iden¬ 
tify  water,  woods,  and  other  (nonwater,  nonwoods)  land  use  (para¬ 
graphs  71-73  and  Tables  12  and  18).  The  unguided  procedure  introduced 

in  this  study  classified  water,  woods,  and  other  land  use  within  5.1  per¬ 
cent  of  the  Satartia  classification  and  within  3.1  percent  of  the  guided 
procedure  for  the  same  data  (Table  18).  This  indicates  that  the  unguided 
procedure  is  adequate  for  most  uses  of  these  three  land-use  classes. 

88.  No  satisfactory  separation  of  inclusive  urban  land  use  was 
obtained  from  the  guided  procedure.  The  results  do  suggest  a  method  of 
locating  possible  urban  areas  (i.e.,  the  red  areas  in  Figures  22  and 

23  identify  surfaces  giving  the  highest  reflectance  values  in  band  4). 
These  are  surfaces  such  as  concrete,  sand,  bare  earth,  rock,  rooftops, 
and  metal,  surfaces  most  often  found  in  urban  areas  and  in  highways  con¬ 
necting  urban  centers.  Though  they  do  not  outline  the  urban  area,  there 
are  aggregates  of  these  pixels  that  suggest  the  locations  of  urban  cen¬ 
ters.  The  results  do  indicate  that  industrial  areas  and  large  business 
areas  can  be  identified  very  well  (Figures  12,  14,  and  33). 

89.  One  of  the  problems  to  be  overcome  in  identifying  a  total 
urban  area  is  the  integration  of  reflectance  from  concrete,  rooftops, 
and  trees  that  make  up  residential  areas  of  the  cities  (paragraph  76). 

90.  This  study  again  attests  to  the  capability  of  Landsat  data  to 
detect  change.  At  least  two  significant  changes  in  land  use  from  1972  to 
1977  were  detected:  the  change  in  an  interstate  highway  through  a  wooded 
patch  and  the  change  in  a  new  residential  area  over  a  5-year  period; 
i.e.,  the  growth  of  trees,  lawns,  and  other  vegetation  dilutes,  to  some 
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degree,  the  high  reflectance  from  rooftops,  streets,  and  driveways 
(paragraphs  82  and  83). 

91.  It  is  evident  that  percent  area  of  the  site  is  not  a  complete 
evaluation  of  the  success  of  the  classification.  The  position  of  the 
pixel  must  also  be  accurate  for  correct  classification  (Tables  18-19). 

92.  The  new  automated  program  for  evaluating  the  classification 
(for  area  and  position)  produced  results  for  area,  as  percent  of  site, 
very  close  to  those  produced  by  using  the  overlay  procedure  (Tables  12 
and  18) . 

93.  There  are  limitations  in  the  use  of  aerial  photography  as 
standards  for  comparison  in  the  success  of  the  Landsat  classification. 
Some  of  these  are  the  availability  of  good  quality  photography,  the 
errors  introduced  in  manual  interpretation  and  digitizing  of  the  land- 
use  patches,  and  the  gridding  procedure  that  may  introduce  differences 
of  a  pixel  width  in  patch  boundaries.  This  has  not  been  researched,  but 
should  be  investigated. 

94.  There  are  also  limitations  in  the  Landsat  data  that  affect 
correct  classification.  Three  important  ones  are  shape,  size,  and  orien¬ 
tation.  The  shape  of  a  land-use  patch  with  very  long  borders  compared 

to  the  total  patch  area  such  as  long,  narrow  patches  increases  the  op¬ 
portunity  for  misclassif ications  when  the  scan  pixel  falls  across  the 
border  and  not  completely  within  or  outside  the  patch.  In  the  large 
patch  with  the  misclassif ication  opportunities  along  borders  (which  have 
an  area  relatively  small  compared  to  the  total  patch  area),  a  much  higher 
percent  of  successful  classifications  is  possible.  The  larger  the  patch, 
the  better  opportunity  for  good  classifications  (paragraphs  80  and  81). 

95.  In  summary,  the  results  of  this  study  show  that  the  guided 
procedure  developed  by  WES  can  be  used  in  sites  other  than  Satartia 
(the  development  site)  to  identify  water,  forested,  and  other  land 
use.  Heavy  industrial  and  commercial  areas  can  be  identified,  but 
the  extent  of  urban  areas  including  residential  parks  and  open  spaces 
cannot  be  identified  with  the  present  system.  Changes  in  land  use  are 
successfully  detected  by  the  classification  system.  An  unguided  pro¬ 
cedure  was  demonstrated  that  proved  to  be  adequate  for  the  three  land-use 
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classifications  in  this  study.  Finally,  automated  procedures  have  been 
developed  to  extract  portions  of  a  Landsat  scene,  rectify  it  to  UTM 
coordinates,  and  fit  it  to  a  gridded  map  derived  from  aerial  photography. 
An  automated  comparison  program  calculates  the  total  area  of  the  maps, 
the  class  areas  in  the  maps,  and  the  percent  of  area  correctly  classified 
in  each  land-use  patch  and  in  the  map  as  a  whole.  The  standard  maps  in 
this  study  were  from  aerial  photography.  The  user  may  select  any  map  of 
his  choice  as  a  standard  for  use  in  the  automated  comparison  program. 


Recommendations 


96.  Based  on  the  findings  of  this  study,  the  following  recommenda¬ 
tions  are  considered  warranted: 

a.  The  effect  of  shape,  size,  and  orientation  of  land-use 
patches  on  Landsat  data  should  be  studied  to  determine 
the  degree  of  accuracy  that  can  be  expected  from  Landsat 
land-use  classification  when  small,  narrow,  or  disori¬ 
ented  patches  exist. 

b.  The  WES  classification  system  should  be  tested  in  sites 
unlike  Satartia  and  Baton  Rouge. 

c.  The  study  of  reflectance  from  urban  areas  should  be 
continued,  perhaps  investigating  variance,  a  definition 
of  texture,  pattern  recognition,  and/or  a  mathematical 
relationship  between  reflectance  bands  of  samples  from 
the  urban  data. 
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Table  8 

Reflectance  Values  From  an  Array  of  Pixels  in  Baton  Rouge 


Band  _ Reflectance  Values ,  CCT  Units 
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Table  10 

Signatures  of  Means  by  Reflectance  Similarity  in  Bands, 
Baton  Rouge  Site,  28  February  1977  (Partial  Listing) 


Band  No.  -  Reflectance 

Value,  CCT  Units 

Count 
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Table  11 

Frequency  of  Occurrence  of  Means  Signatures,  Baton  Rouge  Area 
Landsat  2768-15410,  28  February  1977  (Partial  Listing) 


Band  No.  -  Reflectance 

Values,  CCT  Units 

Count 
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Table  13 

Control  Point  Locations  for  Landsat  and  Aerial  Photograph: 


1 

(Baton  Rouge  Map  Arrays,  1977  Data) 


Control 

Point  Locations 

Standard 

Other 

No. 

X 

y 

X 

Y 

1 

6.30 

8.20 

12.00 

11.00 

2 

14.00 

27.00 

21.00 

29.00 

3 

84.00 

20.00 

90.00 

18.00 

4 

74.00 

63.00 

82.00 

61.00 

5 

9.60 

102.20 

18.00 

107.00 

6 

25.40 

122.00 

36.00 

125.00 

7 

40.70 

81.10 

50.00 

82.00 

8 

90.00 

87.00 

99.00 

85.00 

9 

132.80 

17.70 

137.00 

12.00 

10 

120.70 

57.80 

127.00 

53.00 

11 

104.00 

64.20 

110.00 

62.00 

12 

101.20 

95.20 

108.00 

92.00 

13 

144.80 

66.50 

151-00 

61.00 

14 

125.00 

118.30 

134.00 

114.00 

15 

110.00 

124.60 

119.00 

121.00 

Table  14 

Di fferences  Between  Projected  and  Actual  Distances  of  the  Control  Points  (Grids)  in  Table  13 


co 

r-  £N 

VO 

rH  VO 

vO 

in 

o 

ON 

vO 

cn 

o 

VO 

m 

O  rH  ©  CO  CN 
I  I 


m'v£>oo>r—vj‘oop-<<NCM 


fH  O 

cn 

i 

cn  vo 

i 

ON 

1  1 

m 

ON 

NO 

• 

•  * 

* 

* 

*  * 

* 

■ft 

■ft 

* 

■ft 

* 

* 

o 

1 

© 

1  »n 

CM 

*  * 

* 

* 

*  * 

* 

■ft 

* 

■ft 

* 

■ft 

* 

* 

u  o  w 
D  u 
J3  C 

a  *h 
o  o  c 
u  a, 


Table  15 

Correctly  Identified  Land-Use  Grid  Spaces  for  Baton  Rouge 
1977  Landsat  Classification* 

_ _ Number  of  Occurrences 


Patch  Value,  _ Patch  Value,  Standard  Map 


Other  Map 

1 

2 

3 

4 

TOTAL 

1 

1751 

116 

48 

20 

1,935 

Water 

2 

94 

1862 

470 

148 

2,574 

Woods 

3 

404 

852 

4300 

4143 

9,699 

Open 

4 

58 

71 

166 

1811 

2,106 

Urban 

Total 

2307 

2901 

4984 

6122 

16,314 

Water:  Class  1 

Woods:  Class  2 

Open  space:  Class  3 
Urban:  Class  4 


*  See  paragraph  64  in  the  text. 


Table  16 

Distribution  of  Classes  on  Other  Map  by  Patch  Number  (Partial  Listing) 


Grids  on  Other  Map  Overlaying 


Standard  Map 

Patch 

Patches  on 

Standard 

Map 

Correct 

Class 

Number  of 

Grids  by 

Class 

Number 

Class 

Grids 

1 

3 

4 

5 

1 

5 

5849 

5 

48 

66 

2887 

2848 

2 

1 

1696 

1 

1578 

16 

78 

24 

5 

3 

642 

3 

31 

583 

28 

0 

6 

1 

2 

1 

0 

2 

0 

0 

7 

4 

5 

4 

0 

2 

3 

0 

8 

4 

3 

4 

0 

3 

0 

0 

9 

3 

3 

3 

0 

3 

0 

0 

11 

3 

3 

3 

0 

0 

3 

0 

12 

3 

61 

3 

4 

41 

16 

0 

13 

1 

2 

1 

0 

1 

1 

0 

14 

4 

4445 

4 

38 

58 

4029 

320 

15 

5 

6 

5 

0 

0 

4 

2 

16 

5 

4 

5 

0 

0 

4 

0 

17 

5 

2 

5 

0 

0 

2 

0 

18 

5 

14 

5 

0 

0 

13 

1 

Table  17 


Table  18 


Table  18  (Concluded) 


X? 

U 

<d 

TJ 

d 

<d 

*J 

GO 

ia 

C 

c 

o 

•rH 

4J 

<0 

4J 

0) 

M 

e4 

OJ 


♦o 

<0 

H 


c 

o 

tn 

•H 

U 

If 

u 


<a 

o 

PU 

X> 

c 

03 

a* 

s< 

aj 

u 

u 

< 


o 


u  OJ 
d  4J 
43  -H 
U>  03 
V* 

CD 

a* 


(A 

<U 

M  (A 
5  tA 

o  d 

rH 

o 


tA 

43  43 
43  tA 
U  CA 

43  <0 

H  rH 

CJ 


CA 

M  43 

a;  tA 
.c  ca 
m  d 
O  *— 4 
CJ 


co  cm  cm 


vO 

r- 


on  o 


\p  r*- 

ir> 


xO  xp 


o  co 


fv  N  <N  fO 


r-4  ^  VO  r-i 
ON  00 


co  co  m  i 

r-  in 


in  N  vO  OJ 


h  co  in  n 

ON  00 


00  ^ 

co  oo 


ON  in 
CO  00 


U  rH 

(A 

M  CJ 

<u 

43 

ON 

00 

00 

m 

0) 

43 

CA 

1 

• 

• 

» 

1 

• 

• 

« 

04 

M 

(A 

1 

ON 

m 

1 

00 

00 

NO 

43 

«J 

00 

in 

ON 

r- 

m 

ON 

H 

r-l 

CJ 

-o 

-o 

•H 

•H 

M 

U 

O 

a 

(A 

r- 

in 

rH 

00 

M 

03 

!  E 

• 

• 

E 

j 

• 

• 

43 

(A 

i 

CM 

CM 

l 

o 

CM 

^3 

CA 

o 

CM 

CM 

O 

CO 

in 

/~s 

■U 

CO 

o 

xO 

rH 

in 

VO 

rH 

U  0) 

o 

rH 

r— 1 

r 

•> 

co  ac 
<f)  (0 

CJ 

» 

in 

i 

CM 

m 

x?  a). 

cm! 

C  Ml 

r-* 

CO  o 

CA 

r* 

<r 

rH 

o 

CO 

ON 

NO 

<r 

>4  -<1 

43 

43' 

ON' 

« 

♦ 

« 

♦ 

rH  | 

• 

• 

• 

wi 

43 

CA 

rH 

00 

00 

nO 

>3* 

CM 

VO 

m 

U 

tA 

CO 

r- 

m 

O 

1 

CM 

CO 

o 

co 

J3 

CO 

1  1 

CM 

CO 

m 

CO 

CM 

CO 

m 

00 

H 

rH 

•> 

*. 

d 

r 

r 

CJ 

c 

00 

r^. 

43 

00 

43 

rH  1 

rH 

rH 

r-»| 

<| 

< 

44 

g 

CA 

4J 

00 

ON 

U 

ON 

Q 

M 

43 

M 

• 

• 

o 

• 

> 

43 

tA 

O 

o 

ON 

cv 

ON 

iJ 

43 

CA 

04 

<r 

o 

m 

(Q 

*4 

<0 

vO 

ON 

NO 

o 

4J 

o 

H 

r 

r 

r 

03  y-*v 
U  03, 

CJ 

rH 

m 

rH 

NO 

a  oc 
m  cd 

03  03 
u  M 

c  o 

M  •< 

CA 

CM 

o 

m 

NO 

ON 

rH 

NO 

o 

43 

0) 

• 

• 

• 

• 

CM 

• 

• 

• 

1 1 

43 

CA 

o 

NO 

o 

rH 

NO 

o 

U 

CA 

00 

CO 

CO 

CM 

m 

CM 

vO 

NO 

JC 

si 

<0 

ON 

Os  vO 

ON 

<r 

00 

VO 

a. 

H 

rH 

*» 

•* 

r 

* 

u 

o 0 

r»- 

00 

CA 
0) 
u 
o 

^  c 

<4 

^  M  L  ^ 

Cj  -o  4)  u  a 

H  *)  O  43  P  O 

ie  <0  o  ii 

w  :*  :*  o 

o 

H 


0) 

O 

•rt 

(A 


0) 

u 

co 

cx 

CA 

cs 

a; 


4) 

M 
O 

^  a 

co 

U  CA  W  ,£> 
43  *0  43  U 
H  *J  O  43  P 

ooo** 

J° 

o 

H 


•H 

CA 


43 

U 

«0 

cx 


s 

a* 


(Continued) 


Table  19  (Concluded) 
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STEP  IN  "IDEAL"  PROCEDURE 
STEPS  CURRENTLY  OPERATIONAL 
INFORMATION  FLO* 


GENERALIZED  DIAGRAM  OF 
IDEALIZED  PROGRAMMED  IMAGE 
INTERPRETATION  SYSTEM 


PLATE  1 


DETAILS  OF  SLOPE 
CORRECTION  PROCEDURE 
USED  IN  BLOCK  8,  PLATE  1 


DETERMINE  SLOPE  AND 
ASPECT  ANGLES  OP  THE 
TERRAIN  FACET  CHARAC  - 
TERIZING  EACH  PIXEL 


FROM  EPHEMERIS,  DETERMINE  SOLAR 
ZENITH  AND  AZIMUTH  ANGLES  AT 
TIME  OF  IMAGE  ACQUISITION 


CALCULATE  PIXEL  BY  PIXEL  THE 
SLOPE  ATTENUATION  COEFFICIENT 
.S(*l  -  COS  cul 


WRITE  EQUATION  FOR  LINE 
OF  SIGHT  BETWEEN  PIXEL 
AND  SUN 


COMPARE  ELEVATIONS  OF  ALL 
POINTS  ON  TOPOGRAPHIC  SUR  - 
FACE  (E,l  WITH  ELEVATIONS  OF 
LINE  OF  SIGHT  (E|dsI  TOWARD 
THE  SUN 


DETERMINE  SENSOR  ZENITH 
AND  AZIMUTH  ANGLE 


CALCULATE  cos 


8  10 

WRITE  EQUATION  FOR  LINE 
OF  SIGHT  BETWEEN  PIXEL 
AND  SENSOR 


8  11 

COMPARE  ELEVATIONS  OF  ALL 
POINTS  ON  TOPOGRA  PHIC  SUR  - 
FACE  (E.)  WITH  ELEVATIONS  OF 
LINE  OF  SIGHT  (E,0$l  TOWARD 
THE  SENSOR 


DETAILS  OF  SHADOW 
CORRECTION  PROCEDURE 
USED  IN  BLOCK  8,  PLATE  1 


PLATE  4 


GENERALIZED  FLOW  DIAGRAM  OF 
PROCEDURE  FOR  GUIDED  CLASSIFICATION 


PLATE  5 


REARRANGE  CLUSTER  MATRIX 
TO  GIVE  BEST  APPROXIMATION 
Of  CLEAN  CLUSTERS  (ONE 
TERRAIN  TYPE  PER  CLUSTER) 


SCREEN  CANDIDATE  POINTS 
FOR  NORMALCY 

(5TANOARO  OEVIATION  %  SCREENING 
tINOOVI 


DECISION  00  THE  SELECTEO 
SITES  ADEQUATELY  REPRESENT 
DESIRED  TERRAIN  TYPES’ 


DECISION  IS  IT  LIKELY  THAT 
A  DIFFERENT  SCREENING 
IlNDOf  RILL  SEPARATE 
CLUSTERS’ 


DECISION  IS  IT  LIKELY  THAT  A 
DIFFERENT  SCREENING  UNDO* 

•  ILL  RESULT  IN  A  SNKLE  TERRAIN 
TYPE  IN  EACH  CLUSTER’ 


DECISION  IS  IT  LIKELY  THAT  A  DIFFERENT 
SCREENING  *«00t  «LL  EFFECT  TH| 

OE SIRED  COOOLIOATON  OF  CLUSTERS’ 
(ONE  CLUSTER  PER  TERRAIN  TYPE  I 


t  T 


t  T 


CALCULATE  CCT  CLASS  UMTS  FOR 
EACH  NAVELENGTH  RANO  FDR  EACH 
SELECTED  SAMPLE  SITE  (ROURO  TO 
NEAREST  INTEGER! 


DECISION  IS  IT  LIKELY  THAT 
A  DIFFERENT  VARIANCE  WIN001 
•ILL  SEPARATE  CLUSTERS’ 


DECISION  IS  IT  LIKELY  THAT  A 
DIFFERENT  VARIANCE  fINOOt 
•ILL  RESULT  IN  A  SINGLE  TERRAIN 
TYPE  IN  EACH  CLUSTER’ 


DECISION  8  IT  LIKELY  THAI  A  DIFFERENT 
VARIANCE  flNOOV  IAL  EFFECT  THE 
CONSOLIDATION  OF  CLUSTERS’ 
i ONI  CLUSTER  PER  TERRAIN  TYPE) 


CONSTRUCT  CLUSTER  MATRIX 


DIAGONALIZE  CLUSTER  MATRIX  [ 


DETERMINE  ARE  ALL  {-YESH 
CLUSTERS  DISCRETE’  I 


DETERMINE  DOES  EACH  CLUSTER 
COURT  OF  SITES  REPRESENTING 
A  SINGLE  TERRAIN  TYPE’ 


DETIRMNE  IS  EACH  DESIRED 
TERRAIR  TYPE  REPRESENTED 
BY  A  SINGLE  CLUSTER’ 


DETERMINE  CLUSTER  CLASS  I 
UNITS  IR  EACH  tAVELCWWP 
RAND  | 


8  CLASSIFICATION 
INCLUSIVE  OR  EXCLUSIVE’ 
EXCLUSIVE  '  UNCLASSIFIED  ■ 
PIXELS  ACCEPTABLE 
INCLUSIVE  UNCLASSIFIED 
PIXELS  NOT  ACCEPTABLE 


CLASSIFICATION 
“  INCLUSIVE 


MACHINE  PROCESS.  COOf  ALL 
PIXELS  IN  SCENE  ACCORDING 
TO  STATEO  CRITERIA  (CCT  CLUSTER 
CUSS  LNNTS).  WITH  PROVISION  FOR 
TiNCUSSIFIID  CATEGORY 


DETIRMNE  INCLUSIVE 
CUISTEN  CUB SETS  | 


CODE  AU  PIXELS  IN  SCENE 
(OR  TEST  AREA)  ACCORDING 
TO  APPNOPRIATE  CATEGORY 
OR  TERRAIN  TYPE 


•RITE  OUT  TEST  MAP  OF 
SCENE  (OR  TEST  AREA)  VITH 
EACH  CARDfOATE  CLUSTER 
CLASS  SET 


SELECT  BEST  CLUSTER  CUSS 
SET  ON  BASIS  OF  BIST  NENDITRM 
Of  OVERALL  DISTRIBUTION  OF 
TERRAIN  TYPO 


13  (PLATE  )l 

CLASSIFY  ALL  PIXELS  FOR 
SELECTEO  AREAS  ACCORDING 
TO  CHARACTERIZATION 
SIGNATURE  SETS 


PROCEDURE  FOR  SELECTING 
APPROPRIATE  CLASSIFICATION  CRITERIA 


In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a  facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 


Smith,  Margaret  H 

Acquisition  of  terrain  information  using  Landsat  multi- 
spectral  data;  Report  3:  Application  of  an  interactive 
classification  procedure  in  South  Louisiana  /  by  Margaret 

H.  Smith,  Horton  Struve.  Vicksburg,  Miss.  :  U.  S.  Water¬ 
ways  Experiment  Station  ;  Springfield,  Va.  :  available 
from  National  Technical  Information  Service,  I960. 

87,  [2L]  p.,  [3]  leaves  of  plates  :  ill.  ;  27  cm. 

(Technical  report  -  U.  S.  Army  Engineer  Waterways  Experiment 
Station  ;  M-77-2,  Report  3) 

Prepared  for  Assistant  Secretary  of  the  Army  (RAD), 

Department  of  the  Army,  Washington,  D.  C. ,  under  Project 
LA16i101A91D,  Task  02,  Work  Unit  110-Q6  and  Earth 
Resources  Satellite  Program,  Work  Unit  3158L. 

I.  Data  acquisition.  2.  Landsat  (Satellite).  3.  Multi- 
spectral  data.  L.  Spectrum  analysis.  5-  Terrain 
classification.  I.  Struve,  Horton,  Joint  author.  II.  United 
States.  Assistant  Secretary  of  the  Army  (Research  and 
Development).  III.  Series:  United  States.  Waterways  Experiment 
Station,  Vicksburg,  Miss.  Technical  report  ;  M-77-2,  Report  3. 
TA7.W3L  no. M-77-2  Report  3 


